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ABSTRACT 

In the past several years there has been much Interest in the 
scattering of waves from rough surfaces. Many theories have been 
developed to describe the scattering phenomenon, but these 
generally consider only the return from a single surface and assume 
there is no return from below the surface. This assumption is 
generally violated in nature, which prompted an investigation of the 
effect of a layer on the back scatter signal. 

This work is a study of acoustic wave scattering from a 
layer. The applications of the acoustic results to electromagnetic 
scatter are considered and the limitations are pointed out. 

A theory is developed, using the method of physical optics, 
for tl’ 3 backscattering of acoustic waves from a semi-transparent 
layer having either the front or the rear face randomly rough and the 
other interface plane. It was found that, if the roughness is of the 
order of the incident wavelength or greater, the power returned from 
the two faces could be determined independently and added together 
to obtain the total power return from a layer. 

An extensive experiment was conducted using layers with 
identical roughness and varying thicknesses to compare with the 
theory developed. Each layer was investigated under four 
conditions: (1) rough front and perfectly reflecting smooth rear, 

(2) rough front and smooth rear with both reflection coefficients of 
the same order of magnitude, (3) smooth front and perfectly 
reflecting rough rear, and (4) smooth front and rough rear with both 
reflection coefficients nearly the same. For each case the 
reflection coefficient of the front face was -24.6 dB. Selected 
experimental results were compared with the theoretical scattering 
cross section and found to be in agreement. Comparisons between 
layers of the same configuration and different thicknesses point out 
this increased effect of the rear face as the round-trip attenuation is 
decreased. 

The results of the experiment show that a single-frequency 
system cannot detect the presence of a layer unless the returns from 
the two interfaces can be resolved in range. 
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CHAPTER I 
INTRODUCTION 
1 . 1 Motivation and Purpose 

The problem of scattering of electromagnetic and acoustic waves 
from rough surfaces occurs in remote sensing by radar or sonar systems 
and propagation of waves over rough surfaces. The prediction and esti- 
mation of such scattered signals returned from various terrains is impor- 
tant in designing radar mapping systems and radio altimeters. Radar 
systems in orbit about the earth can yield valuable information about 
vegetation , cultural objects, and geology. For ground-to-ground micro- 
wave communication systems the terrain between the terminals may cause 
scattering of waves . 

Presently available scattering theories are generally based on 
either of the following assumptions* (1) the surface is a perfect reflector, 
(2) in case the rough surface is not a perfect reflector it is homogeneous 
and semi -infinite such that no contribution is received from below the 
surface. While these are very "nice'' mathematically they are quite 
often violated in real life . A somewhat more general approach is to 
consider a lossy layer covering such a surface . Examples of this are: 
a vegetation layer covering the earth , snow covered terrain , soil over 
rock, dust layer on moon or planets, and a sediment layer over a solid 
ocean floor. Experimental evidence of such layers can be found in the 
data presented in papers by Peake (1959) and Taylor (1959). In such 
cases , it seems reasonable to assume the scattering surface to be a 
rough layer covering an imperfectly reflecting surface. Such a layer 
will usually be somewhat lossy. 

Preliminary investigation of backscatter from a layer has been 
conducted by the U . S . Army Engineer Waterways Experiment Station 
(Lundien, 1965) using pulsed radars at P, C, X, and Ka bands. They 
used large laboratory soil samples mounted above flat metal plates; 

. Cl 
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the moisture content, density and thickness of the layer were varied. 

The layer was flat dir both sides and had a maximum thickness of two 

■■'V. 

feet. As tokie expected from such a configuration , the total return 
shows interference effects between the front and rear faces for those 
layers thin enough to have appreciable return from the rear . 

An analysis has also been made for electromagnetic wave scatter 
from a layer with one interface plane and the other rough (Krishen, 1968) 
but only a smooth, perfectly conducting rear face was considered. This 
work also considered acoustic wave scattering with one smooth face 
and one one-dimens ionally rough face. An attempt made to obtain exper- 
imental data for comparison was restricted to the case of a rough front 
only. ' The data obtained were biased, however, by the use of small flat 
facets Which had a slope distribution that peaked away from normal 
incidence . ' 

Four layer models were chosen for this investigation: 

(1) a two-dimens ionally randomly rough interface in front 
and a perfectly reflecting smooth interface in the rear. 


(2) a two-dimensionally Randomly rough interface in front 

and an imperfectly reflecting smooth interface in the rear, ; 

(3) a smooth face in front and a perfectly reflecting two- 
dimensionally randomly rough interface in back, 

(4) a smooth face in front and an imperfectly reflecting two- 

, , dimensionally randomly rough interface in back . 

The scattering from these layers was investigated analytically using the 
Kirchhoff- Huygens method and experimentally using ultrasonic waves in 

Ir 

water. The layers all hacl slight attenuation per unit thickness. The 
results appear in the form of average scattered power densities or average 
differential scattering cross sections . The assumptions used for the 
analysis are discussed in Chapter II. 

The experimental phase of the study used acoustic modeling] 
with carefully controlled parameters . The data available to date on 
layered rough surfaces is meager with surface and subsurface para- 
meters poorly documented. Natural surfaces were usually used and 
measurement is very difficult for the scale of roughness important to the 
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scattering process. The layer used in this experiment is characterized 
by measured thickness and measured statistical description of the 
rough surface. 

1 .2 Discussion of the Kirchhoff Method 

The Kirchhoff method has been used here since the roughness is 
comparable to the incident wavelength . The field scattered by the rough 
interface is formulated according to Huygens' Principle and is given by 
the Helmholtz integral. This integral expresses the field scattered by 
the rough interface in terms of the total field and its normal derivatives, 
or their equivalents , on the surface (Silver, 1941) . The values of these 
two quantities are determined by the tangent plane approximation. The 
criterion for the validity of this approximation is (Brekhovskikh , 1952) 


4 rrpcosQ»^ (ia) 



when the point is not an inflection point. In this inequality p is the 
smaller of the two principle radii of curvature at the point , © is the local 
incident angle and A the incident wavelength. In case the point is an 
inflection point , the condition is . ^ 

Ji 

?4 77 2 cor o » (1.2) 

where x is the coordinate measured along the mean level of the rough 
surface . These requirements restrict this method to locally flat sur- 
faces with small curvatures (gently undulating) . Consequently this 
method cannot be applied in the low frequency limit . 

1.3 Summary of Chapter Development 

In Chapter II an acoustic wave scattering theory is developed 
for a layer with?4ither the front or the back having two-dimensional 
roughness. The assumptions made in the derivation are listed. General 
two-dimensional rough- surface scatter is analyzed with a reflection 
coefficient other than t 1 considered. This theory is then applied to 

■■ ; : ' : ; '■ '' • V.' ■ '' '' ; ''' J 
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the= simplest of the layers considered, the smooth front and the rough 
blck . An expression for the average acoustic power backscattered from 
such a layer is determined . Next the more difficult case of the rough 
front interface is investigated and an expression is developed for the 
average backscatter power from this layer. 


The experimental procedures and equipment are covered in 
Chapter III. The first section describes in some detail the method of 
measuring the average differential Scattering cross section and the 
normalization procedure used in the experiment. A description of the 
acoustic simulation system appears next. This acoustic system was 
used to determine the average backscattered poweiy/bf 1 MHz acoustic 
waves in water ( )\ ■» 1 .5 ) from the layer models used . The target 


was large enough that sufficient independent samples were obtained to 
reduce the averaging err dr to less than 2 dB. The rest of Chapter ill 
describes the acoustic parameters of the material used for the layer 
and the statistical parameters of the rough interface . 

Chapter IV presents the experimental data and a discussion of 
the results obtained. The total average power was obtained for all the 
layers considered and for the thicker layers it was also possible to 
obtain the backscattered power from the two interfaces separately. 

Several comparisons are made between the results obtained for layers 
of the same configuration but differing thickness to illustrate the effect 
of the layer thickness. The effect of the thicker layers is not very 
pronounced, since the signal from the rear face is attenuated. As the 
layer thickness is decreased the effect becomes greater until, for the , 
thinnest, it is not obvious from the return which configuration is being 
observed. Following, the experimental results is a discussion of the 

limitations of the investigation and their effects on the data. The most 
serious limitation is that the smooth surface was not in the far field 
region of the rough surface but any effect of this deviation was not noticed 
in the results. The conclusion reached is that the presence of a layer can- 
not be detected by a single -frequency radar or sonar system unless the 
the range resolution is less than the layer thickness. 

Chapter V contains suggestions for future investigation of the 
nature of scattering from layered media . 



5. 

Background information concerning acoustic waves in solids and 
the effect on reflection from a solid interface is in Appendices A and B 
while a detailed description of the measurements of the acoustic para- 
meters is contained in Appendix C. 

1 .4 The Acoustic-Electromagnetic Analogy 

A) 

The analogy between acoustic waves in water and electromagnetic 
waves in air becomes apparent when the differential equations of each 
case are examined. These "telegraphers equations" appear as (1.3) and 
(1.4). For waves traveling in the plus z-direction they are (Moore, , 



<5 W y _ _ ,, ill 

3 t <j) 2 3 ^ 

where H' is the magnetic field 

E is the electric field 

is the permeability of the medium 

£ is the permittivity of the medium 

p is the acoustic pressure 

U is the particle velocity in the z-direction 
z v 

is the mass density per unit volume 

K is the com pres sibix: ty of the medium 

Combining the two equations in each case as usual yields 
equations of the form 


(1.4) 


wave 
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where c is the velocity of the electromagnetic wave 
em 

c g is the velocity of the acoustic wave in water. 

Thus it is seen that identical equations are valid for the two 
conditions and, as would be expected, the solutions are identical. 
The boundary conditions are 


Electromagnetic 
E tl^ E t2 


Acoustic 

P x = P 2 (1*7) 


H tl - H t2 




( 1 . 8 ) 


where n refers to the normal components and t refers to the tangential 
• ... C) 

components; 1 and 2 refer to the media on opposite Sides of the inter- 
face. The electromagnetic portions of (1.7) and (1.8) show the con- 
tinuity of the tangential components of E and H across the boundary. 
The acoustic condition states that, in an ideal fluid, the pressure on 
both sides of the boundary must be equal and the normal component of 
the velocity must be continuous across the boundary. 

Only the scalar form of the electromagnetic wave equation has 
been considered since an acoustic wave in a fluid is a scalar. The 
restrictions imposed by this limitation have been considered in detail 
by Edison (1960) and Parkins (1965) but basically the analogy holds 

when the vector nature of the electromagnetic wave can be neglected 

= if 

• . If 





il. 


7 . 

and the scalar radar equation is used, and for angles of incidence chat 
are not too large . 
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CHAPTER II 

ACOUSTIC WAVE SCATTERING FROM A LAYER 
2 . 1 Introduction 

The theory of reflection from a one -dimensionally rough layer 
(Krishen and Koepsel, 1968) is here extended to include the scattering 
from a layer with two-dimensional roughness on one face and to allow 
for transmission losses through the layer. The model considered for 
the layer has one smooth and one random, two-dimens ionally rough, 
interface. The model is investigated with first the smooth interface 
and then with the rough one illuminated. The following assumptions 
are made in the derivation of the equations . 

(1) The radius of curvature of the roughness is much greater 
than the wavelength of the incident radiation 

(2) The illumination is a plane wave 

(3) The layer thickness is such that the smooth interface is 
in the reradiation far field region of the rough one 

(4) The observation point is in the far field region of the 
upper interface 

\ 

(5) Shadowing effects on the rough surface are neglected 

(6) Multiple scattering by the rough interface is neglected 

(7) The reflection coefficient is independent of the local 
incident angle over the range of angles of interest. 


2 .2 General Two-Dimensionally Rough Solution 

For a non-perfectly conducting, two-dimensionally rough surface 
state with equation (3.2.1) (Beckmann and Spizzichino , pg. 26) 



(2 . 1 ) 


8 . 
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where 


Figure 


a fa j^s / m - S 't*@ z tc>s ) <K - (5 / *» 0 2 5 1 m ) y 

-fc<?sS>, ■# CoS & z )f] 

k sr (l -*- R ) c<a$$ a -*O~R)c0S 
c = (t + R) $;» &i $*» &* 

JX 

^ 2>* 



/ = f (f } y) is t ^ ie sur f ace height function 
R is the reflection coefficient of the surface 
r is the radius vector y £ +y y -* 
and the angles are defined in Figure 2.1. 


( 2 . 2 ) 


(2 . 3 ) 


(2.4) 


(2.5) 
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If a , b , and c are constant over the illuminated area , i .e . if the 
reflection coefficient does not vary over the range of the local incident 
angle , (2.1) may be integrated by parts to give 



( 2 . 6 ) 


It can be shown that if the illuminated area is much greater than the 
exploring wavelength the only significant contribution comes from the 


first term so (2.6) reduces to 



(2 . 7) 


R f“/ -v cos 6} ces - Six cos <2, 


JL [j 

Cos e, 

L- 


CoS + COS ^2 


( 2 . 8 ) 


A - 4XY is the area projected onto the xy plane 


From (2 , 7) we can compute the average power using the Poynting 

vector 

= zz ( 2 . 9 ) 

i ■ ■ 

where Z is the acoustic impedance at the point and the brackets represent 
the ensemble average of the terms enclosed. The quantity pp* is 



P$_ ffff * t0 ' / y lv~yt ) 
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Defining fiV ti y) 9 f *jLt Z) y t )* nd avera 9 ing gives 




■ m 


i^dr^ + Wyiy-v,.) nJ4, 

e <Ce 


'M 


^ <s/x tf/y W 


/«. 


-v-x-y-v 


( 2 . 11 ) 


Now define the distance between (Xj , and (x 2 , y 2 ) as 
write (2.11) in polar coordinates 


and 


2tT 






Performing the integration over (j) gives 


( 2 . 12 ) 


*> = 




> tr w tr 


(2.13) 


To be considered now is the average of exp(i\^ ) ) which is 

given by 


-oo -«• 

where PC f /t f z ) is the joint probability distribution of the random 
variables ) / and f z . Equation (2 .14) defines the joint characteristic 
function (v tJ -»0 (Davenport and Root, pg. 52) so (2.13) can be 
written as 


Jft 


(2 . 14) 



(2 . 15) 


It is now necessary to assume a specific roughness distribution 
to perform the integration . The most important and probably the most 
typical is the normal distribution . The variance can be chosen to 
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represent any degree of roughness and the correlation distance represents 
the density of irregularities. 

Let / be normally distributed with zero mean and standard 
deviation , in particular let / 1 s distribution be given by 

€*■ 



(2.16) 


Note that this does not say anything about the density of the irregulari- 
ties. This information is conveyed by the correlation function, which 
gives the correlation between the random values assumed by f at 
two points (Xj , y 2 ) and fc 2 / separated by a distance “T . One 
autocorrelation function that can be handled analytically is 

- K ' 

C(tO = e L% (2.17) 


where L is the “correlation distance* for which C('C) drops to the 
_r 

value e 

The two-dimensional normal distribution of two random variables 
j j and ^ 2 with zero mean values, variances <r* and correlated by a 
correlation function C, is 

1 '-f 2 4 2CZ,H t J 2* 


I 






(2.18) 


The characteristic function of this distribution is given by 


X [- v 4 V( i - c)l 


(2 .19) 


Substituting (2.17) into (2 .19) and expanding in an exponential 
series gives 
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( 2 . 20 ) 


For ease of notation let 




or 
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' 2r f (cos + COS <£>,) (2. 21) 


and 

v = Vv 2 ■+ v z 
*y i y , 

(2.22) 


t x 


(2.23) 


Substituting (2.20) into (2.15) yields 

*(**!:■£ 

Z*o L J 

Here the limit of integration has been extended to <*> which is permis- 
sible since for a rough surface the only significant contribution comes 
from the region near f ® 0 . This can also be expressed as 


. L?ttF, 
<£f > - “7T 


09 
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(2.24) 
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by interchanging summation and integration which can be done since 
the convergence of (2 .23) is uniform. From the theory of Bessel functions 


jj'C V r)c mr rdr = £ [R< Cm) > O ) 


(2.25)' 


Applying (2.25) to (2.25) gives 

<2 r" 2 

» L n F 3 -<3 <r 

<4* > - a 5 e M,!M 


*■%* 


(2.26) 


The series converges more rapidly than the exponential series 
but when g is large it is too slow to be practical for computation. To 
enable us to obtain a reasonable solution consider three cases: 

(A) g « 1 > (B) g sssl, and (C) g » l. v 

'a -F 

(A) This corresponds to a slightly rough surface . For this case 


only the first term of the series needs to be considered to give 

.2 


. *. L*n F, q 

< ** > ~ , 


-9 - 




(2.27) 
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(B) This corresponds to a moderately rough surface with the RMS 
roughness comparable to the incident wavelength. In this case (2.26) 
can be used with a reasonable number of terms. 

(C) This is a very rough surface . For this case saddle-point 
integration must be used. Going back to (2.15) and substituting in 
(2.19) and (2.21) we get 

* rwr 




<fr*> = J 1%*) 


(2.28) 


or using (2.17) 

. * 

& 2 TT 

C-f -J 0 > = A* 


yM y *)e tj * 


(2.28)2; 

v- J 


Using only the first two terms of the expansion of exp ^ ?*/£* + S- 

this becomes 


* Srrn Z C ^ ” 

> * — jw?e L TJt, 


(2.30) 


Using (2.25) again we get 
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A* <3 


- «**.* 


e *9 


(2.31) 


or replacing g this becomes 

* _ rr L* F, 1 ' 

AV fi f’ 
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(2.32) 


Thus the avenge scattered acoustic power can be determined for 
any set of 0j , Q ^ • ^ 3 * 
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2.3 Layer with Smooth front Face and Two- Dimensionally Rough 
Back Interface 

This is the simpler of the two cases since the effect of the layer 
can be broken up into a specular signal from the smooth front and a 
scattered signal from the rough back appropriately modified by refraction 
and attenuation . For backscatter the front will not contribute except 
near vertical incidence where it should be the dominant term . For non- 
vertical incidence the main contribution will come from the roughness, 
as modified by the refraction through the front face . 

The incident wave on the rough Surface must be modified to take 
into account the changed incident angle and the signal loss through 
attenuation* This is done by considering the system shown in Figure 2.2. 
Note that © 1 and (ft ^ are related by Snell's law such that 

(f) * s ] * S i u '[« 5 i* 

J (2.33) 


nte* n 


JV--.-. 


Figure 2 . 2 Layer with Smooth Front and Rough Rear 


For the attenuation effect consider a plane wave of unit amplitude inci- 
dent on the top surface . Just below the surface the wave will have the 
amplitude T, i,e . , the incident amplitude multiplied by the transmission 
coefficient. On the rough surface thL'.will then become ' fv 


;<i=T®^ dseo,, i 

;e> j 

The wave scattered back to the sW.wbth surface is given by (2 . 7) as 


(2.34) 
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where (£> ^ has been considered to be zero. Above the layer we have 


(2.35) 


again by Snell's law 
'I = Sh ’ L 


0 - s / m (n s/* 6 > ) 


(2.36) 


and 
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(2.37) 
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Now if we restrict the discussion to the case of backs catter 
which means 




(2.38) 


this becomes 


' F, , 2* 

f‘ ’ i Tr « 


'if 


cJ KcJ 


Y 


(2 . 39) 
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Finding the average backscatter power by the same method as the pre- 
ceding, section yields 
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(2.40) 


where = 


cos 2 6> 


R is the reflection coefficient of the bottom interface 
T is the transmission coefficient from medium 1 to 2 

T” is the transmission coefficient from medium 2 to 1 

g = \J*<r 2 = 4 k* c.os 2 $' o~ z 
= s / ■»' s i* &, ) 
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V r 2/1 ^ 

d = thickness of the layer. 


Rewriting this gives us 

r aTr , - 1 2 - 

<i>f > = 


, O S 


e J * 

l^k^S 'n* 2 $ 
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(2.41) 


which describes the return backscattered from the rear of the layer as 
seen from in front of the layer. 


2 .4 Layer with Twc-Dimensionally Rough Front arjd Plane Rear Face 

This case is much more complex since the rough surface is seen 
twice by the waves returned from the rear of the layer. 

the backscattered power from the rough front is given by (2.26) 
as before with the appropriate choice of parameters and 



R 

COS^ & J 


(2.42) 


To find the field scattered into the layer as a function of direction 
it is necessary to go back to equation (3 . 1 .18a) (Beckmann and Spiz 2 i-* 
chino, pg. 22) and modify it to allow for the transmission. This is 
written as: 

(<? 2 ) - 

where v = k^ - 

p - kj + k 2 ** 

n is the local normal to the surface 
r is the radius vector - xx + yy + f (x,y) *z 





l e 




R, 


(Tt-f)‘Ze JS 


(2.43) 


V 



and the geometry of the system is shown in Figure 2.3. 


% 



Writing the k vectors in terms of the xyz coordinate system 

gives 


kj = ki(sin©i x - cos ©^ 2 ) 


(2.44) 

kg = kg (sin (}>2 cos ^ 3 x + sin 0 g sin $ 3 y - 

cos <p g z) 

(2.45) 

= kg (sin ^g cos 0 3 x + sin^gSin ^ 3 y + cos ^g 2 ) 

(2.46) 

kj = kj (sin 45 g cOs © 3 x + sin@g sin @ 3 

A ' A 

y + cos ©g z 

) (2.47) 

where k^ * u>/v^ and kg » /Vg , 



Using these 


X- 

v - (kjSin^^-kgSin C^cos ^ 3 ) x-kgSin^gS 

in 4> 3 y 

f v 

+ (kgcos ^g-kjCOs© j)z 


(2.48) 
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p = (kjSinOj+kgSin 0 2 cos ^ 3 )x+k 2 sin ^2 sin ^3^ 

- (k 2 cos^ 2 +k 1 cos0 1 )z > (2.49) 

Now consider the factor n db and define the local slope of the 
surface in the x-direction as ft , and in the y-direction as <f . This 
definition gives 

n - sin p cos if x + sin Y cos ft y + cos t cos ft z 
- cos Y cos ft (tan ftx + tan Yy + z) (2 . 50) 

and dS = sec ft sec Y dxdy (2.51) 

such that 

n dS = (tan ft x + tan if y + z) dxdy 

= (JX x + hi $ + z) dxdy '' (2.52) 

3x 3y 

which is the projection of dS onto the xy-^plane. The term (Tv - p) reduces 
to 

Tv - p s[(T-l)k^stn ©j-d+iykgSin ^ 2 cos ^g]x - (T+l)k 2 sin^ 2 sin^ 3 y 
+ [(T+l)k 2 cos </> 2 - (T-l)kjCos e^z 
■s. afc + by + cz (2.53) 


Thus (2.43) can be written as 

i/< g 

i e 

4-n R, 

s ' 

If T is a constant over the range of incident angles encountered during 
any one set of conditions this can be integrated by parts as was done 
for (2 . 1) to give 


e(Q z ) 
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(2.55) 


where the edge effect term is similar to that' ehppuntered before and is 

1 i * * * 

negligable if the illuminated area is large enough . Sblving for the 
bracketed term in front of the integral gives 

c + a ^ + bv y - k a 2 - K, z + T (Kt+KV +2TK fe C os &, * s^0, ^ ^ ) 

^ K 2 cos - k, c©s Q t 




(2.56) 


Thus the field on the rough surface at point Qg due to the scattering 
around is given by 

}} V* 

= l Ar,zX RF - j) e chdy 

-V-Y 

where R Q = d sec . 


(2.57) 


Using this as the incident field in the Helmholtz integral again to find 


iv- r cy/'-r' 



(2.58) 


where / (r -/) #** + ( 1 K,*- 2T'K t K t [c cS e t ' t*< - s » * ^ s; * V c ° * W, J J j 
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T' is the transmission coefficient from medium 1 to medium 2 


To find the average power at we proceed as before to get 


<EO 


= £ 3 rMM e ‘ ^ 1 * *'* <y ~* > * v ‘ u ''' v > * w v ‘' r, \ 
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( 2 . 59 ) 


Making the change of variables 


x - x 

y - y 


x'-x' 


y'-y' 


"T cos ft 
T sin ^ 
T' cos 4 
t'Vin x 


gives 
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( 2 . 60 ) 


Now note that ^ is independent of and (f> of (fi % so this becomes 

<££*> *JCf 

r r' 

| . «> ■ ■■ ( 2 . 61 ) 
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where 
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z + V* 
x T v y 
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Let us now examine the term if Z * * * ' 4 

f; _'•■ 


more detail. This can be written as 


in 




(2.62) 
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where ^ ^ 2 ' ^1 * ' an< * ? 2 ^ are ran ^ om variables with a Gaussian 
distribution function of the form . J . 


W Cz) r 0~i2i t ^ 


*-=- 2<r* 


(2.63) 


where Z is the surface height. 

Equation (2.62) is of the form (Fung 1967) 
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<’ 2 l Z 2 Z 3 Z 4>-' < Z 1 Z 2> < Z 3 Z 4> + < Z lV< Z 2 Z 4> 
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where each of the double averages can be written as 

. r> - *;-■ - • . : f • \ ■' 

i(aS- bS') b z -2*bjs>„.) 

<e > s e 


(2.64) 


(2.65) 


p ss 'being the correlation function between S and S' . Expanding 
the averages in detail gives 
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<r z 1 z 2 > = <- e -i 4 4) > , 

-<r V a 2 C i - fla) 
= e 

(2 . 66a) 

<z 3 z 4 >=< e X'tf.-'iU) 

- o- V c i - ^ 3 .) 

.- e 

(2.66b) 

< Z 1 Z 3 > *<e Ul/i f‘ + l/ * f ‘\ = 

e 

(2.66c) 

< Z 2 Z 4> = 

e -^f-2v;a 4 ) 

(2 . 66d) 

< z i z 4> = <y^ > 3 


(2.66e) 

< Z 2 Z 3> = <e'^^ 4 %r 1 

o 

(2 ..661) 


These are then placed back into (2.61) to yield 
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where 


V/^ sr lK z CoS <f>^ - /<, Cd>S 

vj -.!<i c os ^ -k t cos Q 


6n<i 


vl ' ce>t * 4*z. ~ 2 k t K t c,o%(j> l c*s t )< Z cos Z &, 
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*fV, ' /4 co^ t cos &' + K Z cos*6( 


To find the significance of the correlation functions in (2.67) refer 
back to the geometry in Figure (2.3) and note that the unprimed variables 
refer to the region about Qj and the primed variables refer to the region 
about Qg. In particular, p i2 and ^34 are the surface height auto- 
correlation function about Qj and Q 2 respectively while the rest are the 
correlation of the surface heights between and Q 2 . This means that 
for (J) g near zero (near vertical incidence) the last two terms will 
contribute to the scattered power but that for larger incident angles the 
largest contribution will be from the first part of (2.67) . The first 

part of (2.67) can be integrated as before, (2.13), to give 

a t'i 
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(2.68) 


This then gives the first order scattered power from the layer 
O for large angles of incidence. To better define this note that 

$ 2^3 = 2 d tan ^ 2> 

| For the integrals of the last two terms to be small this distance should 

' be of the order of several times the correlation distance of the surface 

' ■ 0 ■” 

heights . 
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2 . 5 Total Scattered Power from a Layer 

The total scattered power from a layer of either configuration 
considered in this chapter can be determined by simply adding the 
power scattered from the front face to that scattered from the rear face 
of the layer. This is possible in the mean because the scattered signal 

j! 

has random phase. 
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CHAPTER III 

EXPERIMENTAL PROCEDURES 
3.1 Technique of Measurement 

Re-radiation of acoustic waves from a rough surface is described 
by the differential scattehng cross section of the surface. This quan- 
tity is defined by the radar or sonar equation: 


<*«>• 


(3.1) 


where 


is the average of the received power obtained 
over an ensemble of surfaces 

P t is the transmitted power 

is the angle from the antenna bore sight axis 

r. 

G t ( JPJ is the gain function of the transmitter antenna which 
is here assumed to be a function of ^ only 


G rX 


s 


is the gain function of the^ receiving antenna , 
here also assumed to be circularly symmetric 

is the wavelength of the incident radiation 

is the mean surface illuminated 

is the differential scattering cross section. 


The integral of (3 . 1) can be interpreted as follows : the first factor is 
the incident power density, the second factor is the solid angle sub- 
tended by the receiving aperture , and the last factor is the scattering 
cross section. The incident power density is that of a spherical wave 
front at the "surface , modified by the gain-function of the transmitting 
transducer. The scattering cross section was defined in Chapter II 
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on the basis of a plane incident wave front , primarily for mathematical 
simplicity. So long as the correlation distance of the surface heights 
is small enough , the spherical wave front can be considered to be 
locally plane with negligible error. 

Measurement ot the differential scattering cross section can be 
made by using a transmitting or receiving transducer (or both) with a 
narrow beam width . This allows the approximation to be made that the 
scattering cross section is constant over tlie illuminated area and reduces 
the effect of the antenna pattern on the 0 ~ o vs . & curve . For the 
measurement considered here there is no time variation of the incident 
radiation over the illuminated area and the steady state is assumed to 
exist. This is shown in (3.1) by treating P t as a constant and is 
referred to as "beam width limitation . " 

3.! -■ ■-? 

With these restrictions (3 . 1) becomes 


ff SlMM ft ci a 

S (4 ir r*r 


(3.2) 


It is convenient to normalize this average scattered power with respect # 
to the power received upon direct specular reflection from a perfect 
reflector , the smooth air -water interface , when illuminated with P + , 

This power is "J) 

z 

(3.3) 


P e GJo) G r (0) A' 
c ~ (8rrrj) 1 


where t is the distance between the. transducers andr&he surface and 
c 

GJO) and G p (0) are the boresight gains of the antennas. Dividing (3.2) 

• * {/ 

by (3.3) gives 
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where g t ( 4-^ and g r ( •/ J are the normalized antenna gain functions with 
the value of unity for * 1 / » 0 and falling off on either side. The 
integration accounts for the effect of the antenna patterns on the 
scattered signals. 

The determination of <r* 0 for the set of angles is then a measure- 
ment of the average , normalized power and an evaluation of the aperture 
effect for these angles . The antenna pattern functions are given by 
(Morse, 1948) 

2 T,(k a. fi* 4') 

= — 0 . 5 ) 

'Ka 

where Jj is the Bessel function of order one 
k is the wavenumber 
a is the radius of the active face 
.j- is the angle from the boresight axis . 

The transducer used as the transmitter has an active face radius of 8.15 
f0iHmeters which means at 1 MHz, ka - 34 ( A - 1.5 mm) while the 
recfev'er disk corresponds to ka =16. These gain functions are then 
substituted into (3.4) and the integral is evaluated to obtain the aper- 
ture effect . The integration was done by Parkins (1965) using numeri- 
cal methods with the result shown in Figure 3.1. 

The illuminated area must be large enough to permit the surface 
to have a sufficient number of variations to define the random process . 
Clerk (1963) has shown that the dimensions of the illumination should 
be at least 10 times the correlation distance of the surface generating 
^process. For' antennas with 6° half -power beamwidth such as used \ 
here, this requires the rang*: to be greater than 11 times the correlation 
distance, L. 

Also important is the number of independent samples used to 
obtain the average power* Moore (1969) states that the samples can 

(j c 

be considered independent if the value of the correlation between them 
is less than 1/tr , The time corresponding to this decorrelation is - . 

T, = 0 , 65/ fj o fa. 6) 

- . ■■■■■.!. /i* 
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where 

( 5,M °m+<~ Sr " (3 . 7) 

The worst case for this condition is at vertical incidence where , with 
6° antenna beamwidth, 10 mm/sec velocity, and A of 1 .5 mm 

A f d - (.104) = 1.39 Hz, (3.8) 

This means 

T A = .65/1.39 - .467 sec. (3.9) 

Since the averaging time is 20 seconds the number of independent 
samples is 

n t = 20/. 467 = 42.5 (3.10) 

a 

which gives an averaging error of slightly less than 2 dB. 

3.2 Description of the Equipment 

Determination of cr Q requires the measurement of the average 
received power over an ensemble of surfaces . Knowing the antenna 
gain functions for the transducers allows the aperture effect to be 
evaluated. These Measurements were performed using the acoustic 
simulation facility of the Remote Sensing Laboratory of the Center for 
Research, Inc. at the University of Kansas. This facility consists of 
a pair of large water-filled tanks , a mechanical system to provide the 
aiming of the transducers and scanning the surface , and an electronic 
system to generate, receive, and process the ultrasonic signals. An 
overall view of the area is shown in Figure 3 .2 

The transducer aiming system is shown with the transducers 
mounted in Figure 3,3, The angle from the vertical can be changed 
from outside the tank to allow easy determination^ the average power 
vs. angle. Figure 3.4 is a close-up of the transducer assembly show- 
ing the size of the unit and the active areas of the two transducers (the 



Figure 3.2. Overall View of the Acoustic Facility 
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Figure 3.3. Transducer Aiming Mec 'anism 



rigure 3.4. Close-up of Transducer Mount Showing 
Transducers Mounted 



Figure 3.5. Layer in Position with Motion 
Apparatus Shown 




center disks) . The larger transducer is used as the transmitter to 
restrict the illuminated area since it has the narrower beam width, 
while the smaller is used as the receiving unit. 

The method used to scan the surface of the layer is shown in 
Figure 3.5. The target is shown submerged slightly below the surface 
of the water with the transducers below it . The styrofoam blocks on 
the ends of the target are used to support the layer parallel to the 
surface and as a means of connecting the lines used to tow the assembly 
across the illuminated area . The motor and gear box seen on the edge 
of the tank at the upper right of the picture provide the motive force 
to pull the target to the right while a counterweight moves it to the left 
when the motor is reversed. A shunt D.C. motor is used, with field 
and armature voltages adjustable to vary the speed of the scan . For 
each incident angle of the acoustic wave the speed was adjusted to 
allow 20 seconds of averaging to cover the length of the layer. 

The electronic system is shown photographically in Figure 3.6 
and in block diagram form in Figure 3.7. This system is composed of 
four sub-systems: timing, transmitting, receiving, and data processing. 

The timing sub-system is actually two separate units . The first 
part is the pulse generator which synchronizes the entire system , 
establishes the pulse repetition frequency, PRF, and the transmitted 
pulse width. A PRF of 60 Hertz, synchronized to the power line, was 
used for this experiment in order to avoid errors due to the large amount 
of power line hash in the locality . The pulse width was chosen to pro- 
vide a steady-state (beamwidth limited) portion Of the signal from the 
layer at the maximum depression angle . Although it was not essential 
for this experiment the pulse generator is coherently gated by the radio- 
frequency oscillator of the transmitter to reduce the leading edge trans- 
ients of the tone burst . The pulse generator is shown in bljpck diagram 
form in Figure 3.8. The second part of the timing sub-system is the 
delayed sweep portion of a Tektronix 535A oscilloscope which was used 
to select the sample point of the received signal for the processor. 

The screen of the oscilloscope is brightened during the sample gate to 
allow the operator to select the desired time delay. 



Figure 3.6. Photograph of the Electronic System 



Figure 3-7. Block Diagram of Electronic System 



Figure 3-8 . Block diagram of Pulse Generator 
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The transmitting sub-system contains an oscillator which operated 
continuously# an isolating buffer stage, a pulse modulating gate, and 
a commercial wideband, linear power amplifier (Coffman Industries, 

Inc. PA 51). These are all standard circuits and will not be further 
described here. The catput of this sub-system is a radio-frequency tone 
burst of 60 volts peak-to-peak with length determined by the pulse 
generator. 

The receiving sub-system is shown in more detail in Figure 3.9. 
The first amplifier has a high input impedance to avoid loading the \S 
transducer and a 50 ohm output to ‘hatch the step-attenuator following Tt. 
The main amplifier is a commercial unit (Coffman Industries, Inc. , 

PA- 10) with a maximum gain of 40 dB. This amplifier is easily over- 
driven and care must be used to avoid saturation and distortion of the 
waveform . A very high impedance load is needed to avoid loading of 
the output stage . The final amplifier provides this high impedance load 


and boosts the signal with a voltage gain of 10 to reduce the error due 
to detector offset . 


I 


In the data processing system a peak detector rectifies the 
fading RF signals , using the base -emitter junction of a transistor to 
provide a light load on the final amplifier of the receiver. The offset 


voltage introduced by the detector junction is partially cancelled by 
using a common-base , direct coupled amplifier between the detector and 
the sample-hold unit. The sample -hold unit is an Intronics model FS 101 
gated by the +A gate from the Tektronics 535A oscilloscope. The output 
of the sample-hold unit is a series of DC levels — each one the value 


of the detected signal at the end of the sample gate. This output is 
then squared using an Analog Devices model 42 0A wideband multiplier 
to give a series of voltage levels proportional to the power of the 
received signal. These are averaged for a known time interval by an 
operational amplifier connected as an integrator. At the end of the 


averaging period the ^integrator is switched to "hold" automatically and 
the voltage attained is displayed as representing the average power. 
The block diagram of this sub-system is shown in Figure 3.10. The 


processor 


<p> 





(3.11) 
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with T held constant at 20 seconds, "’his calculation of the average 
power was done on a relative basis as the system was normalized with 
respect to a calibration power level. This was the power received by 
the system upon specular reflection from the calm water surface with 
a given transmitted power and a t?tal path length of 1 .5 meters . The 
experimentally recorded power was then normalized as described in 
(3.4) lo eliminate the transducer gains and the transmitter power level 

... ij. * . 

from the calculation of After each set of data the calibration 

procedure was repeated to assure that there had been no gain changes 
during the run. 

3.3 Properties of the Layer Material 

Finding a suitable material from which to form a layer occupied 
considerable time and effort. The properties desired were: ' 

(a) reflection coefficient near 0.5 

(b) low attenuation for longitudinal waves 

(c) high attenuation for shear waves 

(d) able to be formed into the desired shape 

(e) not be adversely affected by submergence ip water 

(f) homogeneous in nature t 

Several materials were tried with none satisfying all the criteria, but a 
soft wax was found which appears to be very close to the desired. The 
major difference with the specifications is that the reflection coefficient 
between the wax and water is quite low. While this reduces the level 
of the return, it is still sufficient to perform the desired experiment. 

One other minor problem is that the wax tends to absorb water after a 
period of time and change characteristics slightly . This can be avoided 
by making all the desired tests within about three days during which 
time no change can be detected . 

3.1,1 Testing the Material 

Preliminary work with the wax consisted of finding a suitable 
mold material to allow the nondestructive removal of the finished target. 
It was discovered early in the testing that the mold should be flexible 
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to allow it to be peeled off and ’'break" the vacuum formed by tho 
close fit. A castable room temperature vulcanizing, rubber proved to 
be satisfactory. To reduce the need. for a mold release agent a silicon 
rubber was used as the mold proper with a backing of lower cost DPR 
#242 rubber to provide the needed support. 

The first acoustic test was made on a slab of wax 37 mm thick 
with both surfaces parallel and smooth. This slab was floated on the 
surface of the water and the acoustic transducers were aimed at it from 
below at normal incidence. This was considered as a very rough test 
to check the reflection coefficient, attenuation and velocity of propaga- 
tion for the longitudinal waves of the wax to determine if it warranted 
further testing. These results were 

R « 0.04 
o Cj « 0.42/cm 
v f ** 2050 m/sec 


v These results coupled with the expectation of high attenuation of the 
shear wave due to viscosity made it a good candidate for further testing. 

The next test was to determine v^ and more accurately. 
To do this a rod of wax was cast with a length of 130 mm and diameter of 
50 mm, enough that the hulk longitudinal velocity was the velocity of 
propagation. With transmitting and receiving transducers mounted on 
opposite ends of this rod the time delay and the received amplitude were 
observed. The same transducers Were then placed in contact with each 
other to obtain the effect of the wax , The procedure yielded: 

! vj =v ; 1870 pi/sec 

oC 9 Q.45 /cm 

J 

which are of the same order of magnitude as the preliminary results . 


•VJ/ 


• . ( 

where R is the reflection coefficient (f 

is the attenuation constant >\ 

r ** v ,’ivi ‘ 

U ' it- 

Vj is the velocity of propagation, 


The shear wave properties remained to be defined after testing 
this rod. To determine the velocity of shear propagation, v„, a prism 
was constructed of the wax. This prism was then submerged in the water 
and illuminated at an angle to excite a shear wave in the wax. This 
shear wave propagated through the wax and emerged as a longitudinal 
wave in the water at an angle determined by the ratio of the velocities. 
Measuring the angle yields this ratio by the prism equation (Born and 
Wolff) . Another feature of this method is that was also determined 
at the same time . This conceptually simple experiment was not very 
successful because a smooth-surfaced prism large enough to occupy all 
the area illuminated by the transducer was not feasible with the limited 
molding facilities available. The results, however, were such that a 
reasonable estimate of vp and v g could be made. These were 

;/v “ 2010 m/sec 
v g = 900 m/sec. 

The results obtained for the longitudinal wave were higher than found 
by the preceding methods but it was expected that the ratio v g /vp 
would be fairly accurate . The main error in the experiment appeared to 
be the curvature of the sides of the prism which would tend to have the 
same effect for both waves. Further investigation using a larger, more 
accurate prism was not undertaken since an easier method of determining 
the velocities was found. 

This method, which was proposed by Hughes (Hughes, et al. , 
1949) required thai a specimen of the material be formed into a circular 
cylinder with longitudinal drive and detection to measure the velocity 
of propagation of both the longitudinal and shear waves. Two such 
cylinders were prepared with the critical dimensions different to provide 

V, |i- 

a check of the results that should be free of any spurious resonance 
effects . The velocities obtained by this method were 
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Cylinder #1 Cylinder #2 

v^ 1800 m/s 1790 m/s 

v_ 500 m/s 552 m/s 

s 

where the agreement between the two cylinders is much better than the 
expected error of measurement. No claim is made to this high accuracy 
but it does tend to add confidence to the method and the results obtained . 

This left only the attenuation constant for the shear wave << s 
to be determined. In order to avoid the inaccuracies that could arise 
from not knowing the mode conversion efficiencies accurately it was 
decided to use a differential method of measurement. To accomplish 
this two different thicknesses of material were illuminated at an angle 
of 66° from normal, beyond the* critical angle for total reflection of the 
longitudinal wave which is 56.5°. This should leave only a shear wave 
to propagate through the slab and be detected. Due to the non-ideal 
antenna pattern of the transducer a small amount of signal was observed 
which had traversed the material as a longitudinal wave, but it was 
sufficiently small that it did not saturate the receiving system to inter- 
fere with the measurement of the shear wave . The resulting difference 
between the strengths shear waves received was assumed to be 
due to the absorption^/ material since all other conditions were 
identical. By this the shear attenuation constant was found to 

be 2.75 neper/cm. 

3.3.2 Results of Tests 


Summing up the results of the various te^ts gives: 


V 

1800 m/s 

v s 

550 m/s 

* v p 

Q ■■ 

9.37 x 10 2 kg/m 3 


, 2.83 x 10 8 n/m 2 


* 2.46xl0 9 n/m 2 
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Here it is to be noted that only the velocities , attenuation 
constants and the density were determined experimentally. The other 
quantities of interest were determined a^s^sfe. in Appendix C using 
the equations developed in Appendi^A % Jffj \ 

The results of thetesfs indicate that, while the wax is not an 

'ideal material from which to form an acoustic layer, it does exhibit 

It 

most of the desirable properties; hence, it was chosen as the material 
to use. k 

f. 

3.3.3 Reflection and Transmission Coefficients vs. Angle 

Using the equations derived in Appendix B and the properties of 
the material given above the reflection and transmission coefficients for 
the cases of water-wax and wax-water interfaces were calculated. 
Figures 3.11, 3.12, and 3.13 show the results for propagation from the 
water onto the wax while figures 3.14, 3.15, and 3,16 are for the other 
case of a wave incident from within the wax. 

3.4 Description of the Rough Surface Used 

The characteristics of the scatter of waves from rough surfaces 
is, in general, determined by the properties of twp random processes. 
These are the variation of surface heights and the variation of the 
material parameters. For most natural or artificial surfaces these pro- 
cesses are of such complexity that they can only be approximated in 
description. This problem Is compounded by i ie lack of statistical 
independence between the processes and their non- statioharity in 
space and/or time , 


.449 

2.65 x 10 9 n/m 2 

8.21 n^r 
.462 n/cm 
2,7 5 n/cm 











51 . 


In order to reduce these problems to a model that can be treated 
analytically, we will deal only with the surface height portion of the 
general case and further restrict it to be a stationary, isotropic process. 
Such a process can be described by its probability distribution function 
of height and the autocorrelation function of the heights as a function 
of separation on the surface. Estimates of these functions were computed 
from a series of sampled height measurements taken along profiles of the 
mold used to produce the rough surface. These measurements were made 
using the profilometer described by Parkins (1965). The probability 
distribution function, autocorrelation function, and the standard deviation 
of heights, <y , were computed from these profiles. Three of such 
profiles were made across separated portions of the mold to allow the 
stationarity of the process to be checked by comparing the sample esti- 
mates from different portions of the surface. 

The estimate of the correlation function was calculated using 
the statistical estimator 



and the sample has a mean value of zero. The denominator of (3.12) is 

— 2 

the sample variance, O' * or the square of the sample standard devia- 
tion. The maximum lag used in the calculation of F was IST/10 as 
recommended by Blackman and Tukey (1958) and N was made sufficiently 
large to adequately define F. According to the sampling theorem (Ben- 
dat, 1958) the samp’ing interval must not exceed half the period of the 
highest frequency component present in the profile. For the surface 
used the value of & x - 0.75 mm amply satisfied this criterion. The 
sample autocorrelation functions computed (using the facilities of the 
Computation Center of the University of Kansas) are shown in Figure 3.17. 





Figure 3.18. Photograph of the Mold Used to Produce 
the Rough Side of the Layers 
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The probability distribution function was determined to have a 
Gaussian nature as expected. The behavior of the sample autocorrelation 
functions shows that the process is quite stationary and smoothly 
undulating as indicated by the close curves and the smooth fall off from 
the origin. It was found that the sample autocorrelation functions are 
closely fit by the Gaussian function 



( 3 . 13 ) 


CHAPTER IV 


EXPERIMENTAL INVESTIGATION OF LAYER SCATTER 

I 

An experiment to determine the effect of a layer on the scattered 
return was conducted using the acoustic simulator described in Section 
3.2. The layer was constructed of the wax with one side being rough 
(see statistical description in section 3 .4) and the other one smooth. 

4.1 Validity of Assumptions Made 

The surface was constructed to fit the assumptions of Chapter II 
as closely as possible. Assumption 2, plane wave illumination, depends 
on the antenna pattern and the range to the surface and was closely fit 
with the transducers used but not exactly. As shown in Figures 3.11 - 
3.16 the reflection and transmission coefficients of the wax-water 
interfacf; change slowly enough with angle to be considered constant 
over the range of local incident angles encountered (assumption 7) . 
Assumption 3 , the two interfaces being in the far field region of each 
other, is violated. To investigate this further a definition of the far 
field region is in order. In the far field region the inequality 


r y 


2j 2 

A 


where R is the distance from the surface to the observation point 
is the length of the scattering area and 
A is the incident wavelength, 

must hold. Barrick (1965) has shown that this condition can be replaced 


by 




where L is the surface correlation distance for a Gaussian correlated 
surface. This condition means the layer thickness must be 

d > y R 

54. .■:> . 
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for the far field assumption to be valid. For the surface used in the 
experiment (L ** 11.6 mm) this becomes 


d 2(11. 6) 2 . 2(135) 
a * 1.5 " 1.5 


- 18.0 cm 


This thickness was such that no return could be obtained from the rear 
face. 


4.2 Theoretical Average Scattering Cross Section 

The theory developed in Chapter II points out that the average 
power from the two interfaces can be determined separately and added 
together to obtain the total average power from a layer. 

The expected return from the smooth surface is that of a specular 
reflector with the proper reflection coefficient. This is a convolution of 
a delta function, the specular reflection of a plane wave, and the 
antenna pattern. The effective antenna power gain was obtained experi- 
mentally by recording the received power from the smooth water-air 
interface as a function of incident angle . This was then normalized to 
unit area using the aperture effect to obtain the result given in Figure 

4.1 • 

The signal from the rough interface is described in terms of the 
variation of average scattering cross section with incident angle. 
Equeition (2.26) was solved using the measured parameters of the surface 
roughness to obtain the curve presented in Figure 4.2. This, when 
scaled as above, describes the contribution of the front face for the 
rough first case directly. For the smooth first case this must be modi- 
fied to take into account refraction at the smooth interface as was done 
for Figure 4.3. 

These curves can be combined by adding together, with the 
appropriate scaling for attenuation, reflection and transmission coeffi- 
cients, to obtain the total average scattering cross section of the layer. 

4.3 Experimental Results 

The four configurations of the layer mentioned in Chapter I were 
examined in detail as a function of layer thickness and attenuation. 
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Theoretical Scattering Cross Section of the Rough Surface as 
Modified to Account for Diffraction Through the Smooth Front. 
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All the layers were prepared from the same mold in order to have identi- 
cal roughness as the thickness was varied. The results of individual 
layers are presented first then a comparison of three thicknesses is 
offered for each configuration. Following these is a comparison of results 
from various configurations. 

4.3.1 Rough Front Face and Perfectly Reflecting Smooth Second Face 

For this configuration the Fresnel reflection coefficient for the 

-2 

front face is 5.9 x 10 at normal incidence and for the tear face is 
unity for all angles. These were obtained by floating wax on the 
surface of the water such that the first interface was water-wax and the 
second was wax-air. The layer thickness used for this condition ranged 
from 11 cm to 1 cm. The thickest layer. Figure 4.4, did not exhibit any 
layer effect due to the large amount of total attenuation of the signal 
from the rear face. As the thickness of the layer was reduced the specu- 
lar term from the smooth rear face becomes more pronounced until for 
the 1 cm layer. Figure 4 . 14 , it was the dominant return . The complete 
set of curves for this configuration is presented in Figures 4.4 through 
4*14. Figure 4.15 is a composite of the total average scattering cross 
section of the 11 , 4.5 , and 1 cm thick layers to point out the transition 
with thickness . It is noted that the effect of the rear face is evident 
only near vertical incidence implying that the average power represented 
by (2,68) is negligible compared to the rest of the terms. 


4.3.2 Rough Front Face with Imperfectly Reflecting Smooth Rear 

As above the Fresnel reflection coefficient of the front face was 

-2 r -2 
5.9 x 10 but now the rear face was 5.7 x 10 . This was accomplished 

by weighting the wax to sink it below the surface to have water-wax 

then wax-water interfaces . The target was supported parallel to the 

surface of the water by styrofoam blocks . The thickest layer to exhibit 

a layer effect in this configuration was 3.5 cm. Figure 4MIQ shows this 

as a very slight peaking of the curve hear vertical incidence. The rest 

of the layers, Figures 4.17 through 4.21 , exhibit more of this effect 

but never to the same extent as the case with a perfectly reflecting 


rear face. The 1 cm layer for this configuration is quite similar to the 

f 

3 cm layer with a perfectly reflecting, smooth rear face. Other then the 
reduction of the magnitude of the return from the rear face the scattering 
cross section curves are similar to those obtained for the perfectly 

Tr 

reflecting case . A composite curve of the total scattering cross section 
for the layers of thicknesses of 3.5, 2.5, and 1 cm is given in Figure 

i 

4.22. With both sides having nearly the same, low reflection coefficient 
it is noted that the effect of the' layer is slight even for a. thin layer. 


4.3.3 Smooth Front Face and Perfectly Reflecting Rough Rear Face 

The layers used for this configuration varied from 11 cm to 1 cm 

as before . The normal incidence reflection coefficient for the front 

-2 <• 
face was 5.9 x 10 and for the rear face it was unity for all angles. No 

- '-cycr effect was noted for the thickest layer. Figure 4.23 and the fCT - 

- .. - • j- P 

■v-vr - 0 curve is very similar to that obtained from the water-air interface. 
There is evidence of small scale roughness on the bas i caliy flat surface 

\ _ r ' * 

as the drop off from vertical incidence is not quite so’ rapid as the* antenna 

\ * 

power curve. As the layer thickness was reduced'the effect of the backr 
surface roughness becan to appear -- first asu "further lifting of tfte S" 
curves for the large incident angles and finally, for a thickness of .1 
4.5 cm (Figure 4.2 7), even at normal incidence. For the layers thinner ; 
than 2 cm , Figure 4.31 , the return from the front face is completely 
masked by the effect of the highly reflecting rear face. This transition . • 
is quite evident in Figure 4*34 which compares the cT 0 vs . & curves 
for the 11, 4.5, and 1 cm layers . 

4.3.4 Smooth 'front Face and Imperfectly Reflecting Rough Rear Face 

Here the layer was again submerged to give a normal incidence 

-2 ' *. 

reflection coefficient of 5 . 7 x 10 from the rough rear face. A ret-urn 
from the rear face was -first noted for the 4.5 cm thick layer where lit had 
no effect near the vertical but became important for incident angles 
greater than 15° (figure 4.35). This pattern continued as the thickness.; 
was reduced .until for the 2 cmj.ayer (Figure 4 ._3j9) some. effect is noted 
at. vertical incidence.. Further reduction of the ’thickness' continued 
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Figure 4-4. Scattering Cross Section. of Layer with Front Rough , R-.059, 
and Rear Smooth, p = i . 0 , for Round-Trip Attenuation of 88 dB. 
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Figure 4-5 . Scattering Cross Section of Layer with Rough Front, R -.059 , 
and Smooth Rear, R=1 .0, for Round-Trip Attenuation of 72 dB. 
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Figure 4-7. Scattering Cross Section of Layer with Rough Front, R=.059 , 
and Smooth Rear , R=1> 0, for Round-Trip Attenuation of 44 dB. 
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Figure 4-8. Scattering Cross Section of Layer with Rough Front, R=,059, 
and Smooth Rear, R-1.0, for Round-Trip Attenuation of 36 dB. 


SCATTERIN 



INCIDENT ANGLE 

Figure 4-9. Scattering Cross Section of Layer with Rough Front, R=.059, 
and Smooth Rear, R=1 .0 , for Round-Trip Attenuation of 30 dB. 
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Figure 4-12. Scattering Cross Section of Layer with Rough Front, R = .059, 
and Smooth Rear, R=l. 0, for Round-Trip Attenuation of 16dB. 




and Smooth Rear, R=1 .0, for Round-Trip Attenuation of !3dB. 
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Figure 4-14. Scattering Cross Section of Layer with Rough Front, R-.059, 
and Smooth Rear, R=1 .0 , for Round-Trip Attenuation of 8 dB. 
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Figure 4-16. Scattering Cross Section of Layer 

and Smooth Rear, R=. 057, for Round- 


Rough Front, R=. 

Trip Attenuation of 30 dB 
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Figure 4-18. Scattering Cross Section of Layer with Rough Front, R=.059, 
and Smooth Rear, R-. 057 , for Round-Trip Attenuation of 2 0 dB* 
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Figure 4-19. Scattering Cross Section of Layer with Rough Front, R=r.059, 
and Smooth Rear, R-.057, for Round-Trip Attenuation of !6dB. 
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Figure 4-20. Scattering Cross Section of Layer with Rough Front,, R=.0S9, 
and Smooth Rear, R = .057, for Round-Trip Attenuation of 12dB. 
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Figure 4-22. Scattering Cross Section of Layers with Rough Front, R^.059, 
and Smooth Rear, R-. 057, as a Function of Attenuation. 
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and Rough Rear, R=1 .0, for Round-Trip Attenuation of 88 dB. 
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Figure 4-25. Scattering Cross Section of Layer with Smooth Front, R=.0S9 
and Rough Rear, R-1.0, for Round-Trip Attenuation of 56 dB. 
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Figure 4-26. Scattering Cross Section of Layer with Smooth Front, R-.059, 
and Rough Rear, R-l .0, for Round-Trip Attenuation of 44 dB. 
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Figure 4-27 . Scattering Cross Section of Layer with Smooth Front, R = . 059 , 
and Raugh Rear, R=1.0, for Round-Trip Attenuation of 36 dB. 
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Figure 4-28. Scattering Cross Section of Layer with Smooth Front, R=.0S9 
and Rough Rear, R-1.0, for Round-Trip Attenuation of 30 dB. 
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Figure 4-29. Scattering Cross Section of Layer with Smooth Front, R=.059, 
and Rough Rear, R»1 .0, for Round-Trip Attenuation of 24 dB. 
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Figure 4-30. Scattering Cross Section of Layer with Smooth Front, R=.059, 
and Rough Rear, R=1 .0, for Round-Trip Attenuation of 20 dB. 
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Figure 4-32. Scattering Cross Section of Layer with Smooth Front, R=.0S9, 
and Rough Rear, R-1.0, for Round-Trip Attenuation of 12dB. 
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Figure 4-33. Scattering Cross Section of Layer with Smooth Front, R-.059, 
and Rough Rear, R-1.0, for Round-Trip Attenuation of 8 dB. 
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Figure 4-34. Scattering Cross Section of layers with Smooth Front, 
and Rough Rear, R-l .0, as a Function of Attenuation. 
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Figure 4-35. Scattering Cross Section of Layer with Smooth Front, R=.059, 
and Rough Rear, R = .057, for Round-Trip Attenuation of 36 dB. 
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Figure 4-37. Scattering Cross Section of Layer with Smooth Front, R-. 059, 
and Rough Rear, R-.057, for Round-Trip Attenuation of 24 dB. 
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Figure 4-38. Scattering Cross Section of Layer with Smooth Front, R=.059, 
and Rough Rear, R=.057, for Round-Trip Attenuation of 20 dB. 
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Figure 4-39. Scattering Cross Section of Layer with Smooth Front, R=.059, 
and Rough Rear, R-.057, for Round-Trip Attenuation of 16dB. 
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Figure 4-40. Scattering Cross Section of Layer with Smooth Front, R=.059, 
and Rough Rear, R=.057, for Round-Trip Attenuation of 12dB. 
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Figure 4-41. Scattering Cross Section of Layer with Smooth Front, R=.059, 
and Rough Rear, R = . 057 , for Round-Trip Attenuation of 8 dB. 
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Figure 4-42. Scattering Cross Section of Layers with Smooth Front, Rr.059, 
and Rough Rear, R-.057, as a Function of Attenuation . 
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the trend toward a return resembling a rough surface with a peak near 
the vertical. Figure 4.42 shows this trend as the layer thickness is 
changed. 

4.4 Comparison of Experimental and Theoretical Results 

Since this is primarily an experimental investigation no detailed 
comparison of the experiment and the theory will be attempted . The 
simpler cases can be compared directly with the three theoretical 
curves , Figures 4 . 1 through 4.3, however. 

4.4.1 Scattering Cross Section of the Smooth Face Compared to 
the Antenna Power Gain 

For the thickest layer with a smooth front face the return from 
the rough rear was attenuated enough that it had no observable effect 
on the return. This case was used to examine the scattering cross 
section of the smooth side of the layer by superimposing Figure 4.23 
and Figure 4 .1 , the antenna power gain curve , in Figure 4 . 43 . The 
experimental curve fits the antenna curve very closely from vertical 
down to 6° . For larger incident angles the return from the smooth side 
of the layer is much more than expected, indicating the presence of 
small scale roughness on the supposedly plane surface. Visual obser- 
vation of the flat side of the layer revealed just such an effect in the 
form of a slight "aliigatoring" of the surface. 

4.4.2 Scattering Cross Section of the Rough Face Compared to the 
Theory 

Again the thickest layer did not have any return from the rear 
face when the rough side was illuminated first so it was used. Figure 
4.44 is a plot of the theoretical scattering cross section with the experi- 
mentally determined points shown for comparison. There is a deviation 
of 2 dB at vertical incidence but the rest of the curve is fit quite closely. 

4.4.3 Scattering Cross Section of the Rough Face in the Rear 
Compared to the Theory 

The thinnest layer with the rough perfectly reflecting rear face 
was used for the experimental data since the expected effect of the 
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Figure 4-43 . Comparison of the Scattering Cross Section of the Smooth Side 
and the Antenna Power Gain Function. 
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Figure 4-44. Comparison of Experimental and Theoretical Scattering Cross 
Sections of the Rough Surface Used. 
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smooth front is 16 dB below the return from the rear at vertical incidence 
and even lower as the angle is increased. As seen in Figure 4.45, 
the experimental points are very near that predicted by the theory, 
except for vertical incidence, as expected due to the effect of the smooth 
front. 

4.5 Summary £ t..e Experiment and Conclusion : 

The experimental layer model used had one smooth face and the 
other two-dimensionally rough. Both rough front and smooth front cases 
were considered with the front face having a reflection coefficient of 
0.059 for all configurations. The rear face was either perfectly reflecting 
or had a reflection coefficient of 0.057. For each case the layer thick- 
ness was varied from that for which no return was received from the rear 
to the thinnest that could be produced. The summary of each configura- 
tion as a function of round-trip attenuation in the layer, layer thickness , 
is given in Figure 4.46. For each case it is noted that the effect of 
the rear face • " small for the layers with high total attenuation but 
becomes increasinbly important as the layer attenuation is reduced,. 

Also of interest is the effect of the configurations with the same 
attenuation within the layer. This is shown in Figure 4.47 where four 
sets of curves are presented with the attenuation fixed for each set. 

The most important observation here is that the layer effect is greater 
when the reflection coefficient of the rear face is high than when it is 
low. Also noted is that, when the layer attenuation is low and the 
reflection coefficient of each face is of the same order of magnitude, 
the return is basically the same no matter which face is illuminated 
first (see curves 2 and 4 of Figure 4.47a). 

Another comparison to be made is one in which the layer attenua- 
tion and the reflection coefficient of the rear face are such that the 
return from the rear is expected to be the same for two cases. For the 
rear reflection coefficient of 0.057 (-24.8 dB) and 8 dB layer attenua- 
tion the total effective attenuation is 32 .8 dB. This could be compared 
to the case of a perfectly reflecting rear face with 32.8 dB attenuation 
in the layer. No such layer was examined with the closest being one 
in which the layer attenuation was 30 dB, which should still be a good 
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(c) Attenuation = 3 6 dB 

1 Front Smooth R-.059 
Rear Rough R - 1 . 0 

2 Front Smooth R=.059 
Rear Rough R = .057 



(d) Attenuation = 88 dB 

3 Front Rough R-.059 

Rear Smooth R'1.0 

4 Front Rough R-.059 

Rear Smooth R-.057 


Figure 4-47. Comparison of the Four Configurations with Fixed Attenuation 
in the Layer . 
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Figure 4-48. Comparison of TwO Layers with Smooth Front Face and the 
Same Expected Effect of the Rough Rear Face. 
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comparison. When these two thicknesses are compared for the layer 
with the .smooth side first, Figure 4.48, it is seen that there is very 
little difference in the curves except for incident angles beyond 15° 
where they start to separate. Even for the extreme of 35° incident angle 
the difference is only 2.5 dB which is within the total error possible in 
the measurements . 

The effects, if any, of the violation of the assumption that the 
smooth face is in the far field region of the rough face should be more 
noticeable when the rough face is illuminated first. Making the same 
comparison as above except illuminating the rough face first. Figure 
4.49 , fails to show any change in the shape of the curves . Thus, for 
the roughness and materials used in this experiment, there appears to 
be no effect due changing the layer thickness except that due to the 
attenuation within the layer. 

The most important goal of this investigation was to determine if 
the presence of a layer could be detected by backscatter investigation of 
an unknown surface. The conclusion is that it cannot be done with a 
single frequency system unless the range resolution is fine enough to 
allow the separation of the return from the two faces. This can be done 
using a system with a very narrow antenna pattern at vertical incidence 
and a range resolution smaller than the layer thickness . 

The implications of this conclusion are quite important. It 
casts doubt on the validity of the curve fitting procedures used to pre- 
dict the roughness of terrain, both terrestrial and celestial, from radar 
backscatter measurements. The majority of these measurements are 
made with no knowledge of the presence or absence of a layer. Refer- 
ring to Figures 4.46 and 4.47 one can see that even for the same surface 
conditions the return can be greatly affected by a layer. However, if 
addionai information is available, the presence and nature of a layer may 
be remotely sensed. For example, comparing returns from the "ame 
terrain in winter and summer could yield the depth and moisture content 
of a snow blanket. 
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CHAPTER V 


RECOMMENDATIONS FOR FUTURE INVESTIGATIONS 


This work was undertaken to provide a starting point for the study 
of a general rough layer. As such it was restricted to a single frequency 
investigation of a particular layer situation. A logical extension would 
be to repeat the same experiment for different frequencies to see if the 
presence of a layer can be detected in this manner. After this the layer 
material and roughness parameters should be varied systematically to 
determine their effect . 

Extending the theory to include a layer with both surfaces rough 
is feasible even though it will generate complicated multiple .integrals 
which may not be possible to evaluate analytically. The experimental 
investigation of a layer with two rough surfaces should be conducted 


The effect of inhomogeneities within the layer may be of equal, 
or even greater, importance for practical cases; hence, it should be 
investigated thoroughly. 

Polarization effects may be quite important. For example, the 
Brewster angle for front and back faces (present only for “vertical" 
polarization) may be different arid cause serious changes in the shape 
of the curves . Ultrasonic simulation , of course , cannot show these 
phenomena, so electromagnetic experiments are called for. 
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SOUND PROPAGATION IN AN EXTENDED ELASTIC MEDIUM 


A-l . Stress and Strain on an Element 


The analysis of this appendix follows that of Kolsky (1953) with 

part of the notations changed to be more consistent with the usage currently 

being employed in acoustic simulation. 

The stress on a surface element in a solid body has components 

both normal anil /tangential to the plane. Using three mutually orthogonal 

axes, O , O . O , nine components of stress act on three planes normal 
x y z I 

to these axes and passing through a f'oint p. These will be denoted by 


^xy ‘ 


yy ' w 22 ' ^xz ' <r xy 
script denotes the direction of thit otr^'s and the second letter defines the 

plane in which 'It acts. By considering an infinitesimal rectangular para- 

llelopiped around P with its faces normal to the axes, as shown in 

Figure 1, and taking moments , it is seen that for static equilibrium 


<r xy 


The first letter in the sub- 


crjy- &~yx' ^xz- ^zx’ ^yz zy' so that only six independent 
components of stress remain. The stresses acting on any other plane 
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The displacement of any point V in the body may be resolved into 
components u, v, and w along the x, y , and z axes such that if the ori- 
ginal coordinates of P were (x, y , z) they would become (x + u, y + v, 
z + w) . The strain at the point is defined by considering how its position 
relative to adjacent points has changed. Consider a point very close to 
P, with undisturbed coordinates (x + 5 x, y + 6 y , z + 6 z) and let the 
displacement it has undergone by(u+ 6 u, v + 5v,w + 6w) . These can 
be expanded by a Taylor expansion of the form: r Q, 



u * 



l. la) 


/XT 




(A. lb) 




2 W' , 




2 w 
2* 


Si 


* » • 


(A.lc) 


For sufficiently small Sx, S V > and £ 2 only the first order terms are 
significant and these become: 


■ft' 

3 

i* 

. 2* . 
T~*Y 

. 3*£ J** 

(A. 2a) 




3 a 

(A. 2b) 

f*- 

f it* 



(A. 2c) 


Thus if the nine partial derivatives in (A. 2) are known at the point , the 
relative displacement of all surrounding points may be found . These ni^> } 
"Quantities may be regrouped and defined as: if 
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£ - **- 
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dl 
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_ dw 3\/ 


2>y ai 
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The first three quantities £ 


xx' 
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av 

jp 

yy- 


2 u> 


s (A. 3) 


ax 


‘y, 


£ „ are seen to correspond 
z z 


to the fractional expansions and contractions of line elements passing 
through P parallel, respectively, to the x, y, and z axes. The second 
three £ x , € , and € correspond to components of sheer strain in 

the planes denoted by their subscripts. The last three u-> x , tv y # andu> z 
do hot corresp^:Xd to a deformation of the element but are the components 

, *.• - F \y' 7 * ’v - ' \ 

't.- ’ *i* - ‘ J . f- i 

of its rotatl#X%s a rigid body. Note that if x, 'y , and are unit vectors 
in the respective directions 


curl s 1 ( + <L*)y +. ( iit 

V*y >«./ W dir ^di ay 


)* 


so that if the displacement of the element is regarded as a vector, "J, 
the components of j curl sf are o^ , co y , and uj z , 

The first six quantities in (A. 3) are called the components of 
strain . When the last three are zero , curl s = 0 , the deformation is 
irrotational and described as pure strain. 


o 


- /' 


A-2 . Generalized Hooke's Law 

Within the elastic range of the material each of the six components 
of stress at any point is a linear function pf the six-components of strain. 
This may be expressed as: a 


f" ", 




<r- * 
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(A. 4) 


where the coefficients are the elastic constants of the material . 

In order for the elastic energy to be a univalued function of the 

strain any coefficient c must be equal to c_ (Love, 1927). This reduces 

rs sr 

the number of independent coefficients from 36 to 21 . In an isotropic 
solid the values of the coefficients must be independent of the choice of 
axes . This reduces the independent constants to two , namely and • 
Then 

= A 


c 12* c 13* C 21 G 31 ” °23 " c 32 


V*3 


c 44 " C 5S = c 66 * K 


v, C U - c 22 * c 33~Ai + 2 ^»" 

ft 

and the., other 24 coefficients become zero; Equation (A. 4) may thus be 
written as: 


“ A* & 4 
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where 


A=^ + e,v* e » = r * 4 77 * If =d.vi 


The change in volume of a unit cube represented by A is known as the 
cubic dilatation. 

The two elastic constants, A m and / *4 m , are known as Lame's 
constants and completely define the elastic behavior of an isotropic 
solid. Four other elastic constants are commonly used. These are 
Young's modulus, E, Poisson's ratio, V , the bulk modulus, k, and the 
rigidity modulus which is identical withyM m - Using (A. 5), E, y , and k 
may be expressed in terms of and^ m . E is defined as the ratio 
between the uniformly applied plane stress and elongation with the lateral 
surface free from constraint. Consider the case of applied and all 

other stresses being zero. The first three equations then become: 


e ■*-. 
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e yy 

+ e V2 

c +■ 
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Solving for £ xx , € , and £ zz it is seen 
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Young's modulus , E , is given by tr^/ £ x5F; so: 
yU** 0^ ^ *»» + m ) 


E = 




(A. 7) 
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Poisson's ratio, V* , is defined as the ratio between the lateral contraption 

and the longitudinal extension, i ,e. - £ / £ . Hence 

YY 


m 


V s 




f]\ 


(A. 8 ) 


The bulk modulus, k, is defined as the ratio between applied pressure and 
the fractional change in volume when the specimen is subjected to hydro- 
static pressure. Under these conditions jr = /r* 0 -* = Pe 

xx yy zz 

<r yz = ^ xz - <r xy - 0 so that from (A. 5) 


-P 


and the fractional change in volume is 


(A. 9) 


so that 
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A-3 . Equations of Motion 

In order to obtain the equations of motion for an elastic medium 
consider the variation in stress across a small parallelopipedon with the 
. sides parallel to the set of rectangular axes , Figure 2. The components 
of stress will vary linearly across the faces so to obtain the force acting 
on each face the r value of the stress at the center of each face is multi- 
plied by the area of that face . 


tt' 

’A, 

\\ 








Figure 2 . Stresses on Parallelopiped 


As seen from the figure , six separate forces will be acting parallel 
to each axis . If the resultant force in the x-direction is considered it 
becomes: 


Ik* 


6 X 

which simplifies to 


J(<r + IX^Kcr- * 


f^X 

V > v' 
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If the body forces such as gravity are neglected the application of Newton's 
second law of motion shows that this must equal 



l4o 


where p is the density of the element and u is the displacement in the 
x-direction as before. This gives 


Y dt a* £y - 
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(A. 15a) 


and similarly for v and w. 
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These equations of motion hold no matter what the elastic constants 
of the medium are but in order tp solve them for a specific case the elastic 
relations must be used. For ah isotropic medium these are given by (A. 5), 
and substituting from these for the stress components in (A. 15a) gives 


eff* - + §y Uej+ k u,f„ I 


(A. 16) 


Now from the definitions of (A. 3) 


e * — a 


Thus (A. 16) becomes 
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Rearranging terms and using the relation, V s ( ? — + - — ) 
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Similarly 
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Differentiating (A. 19) with respect to x gives 


a*n 


P bt dt z " 


a 4 
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Repeating for (A. 20) and (A. 13) 
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adding these become: 


p if* = K +2 /0^ 2 A. 


(A. 2 3) 


(A. 24) 


This is a wave equation and shows that the dilatation, is 
propagated through the medium with a velocity of ♦ ^>*)/pY z - 
As another choice eliminate between (A. 20) and (A.21). 

P £ If ■ (K^.)& 
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Subtracting 


3* / 3^ _ 3^ 
P 3t z ( dv 3* 
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But = Z a, 

dy d*. x 

a 2 a,. 


by (A. 3) so this can be written as 


0 —1 - aJ XT 2 
K dt‘ MJ x . 


(A. 28) 


Similar equations may be obtained for cu and u> 2 . Thus the rotation is 

-m— pm ■* 

propagated with velocity ^ . 

If the dilation is zero {A. 15), (A. 16) # and (A. 19) become 


n ^ t V7 2 
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The condition that the rotations u^, u^,, and to z all vanish is 
u, v, and w satisfy the conditions 


(A. 29c) 
if 
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where is a potential function. Then 


A ~ V a 0 and 


a a 

d* 


^ V M 


Substituting these into {A. 15): 


p ^5 
o ^ 


(A. 30) 


and similarly for v and w . 
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Thus it is shown that in the interior of an elastic solid waves may 
be propagated with one of two different velocities. Waves involving no 

rotation travel with velocity c^t ■* 2 while waves involving 

no dilatation travel at c 2 = fp. Strictly speaking these waves should 
be called irrotatlonal and eguivoluminal respectively. The irrotational 
waves are more commonly called dilatational , longitudinal, com pres sional, 
or P-waves while the equivolurninal waves are known as transverse, 
lateral, or shear waves . In keeping with most of the literature of sound 
propagation the remainder of this writing will refer to them as the longi- 
tudinal and shear waves. 

It can be shown that any plane wave propagated through an iso- 
tropic elastic medium must travel with one or the other of these velocities , 
c^ or c 2 . Since the medium is considered isotropic there is no loss of 
generality if a plane wave propagating parallel to the x-axis is considered. 
Let its velocity of propagation be c; then the displacements u, v, and w 
will be functions of a single parameter - x - ct. 

’ Then: •'/ 
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and the partials with respect to y and z are all zero. Substituting these 
into the first equation of motion {A. 19) gives: 

- (^vv, V**) fv**. ' 




(A. 32) 


Similarly from (A.20)and (A. 21) 
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(A. 34) 


Equations (A. 32), (A. 33), and (A. 34) can be satisfied in one of two .ways: 
either c 2 =(A m *2 /J .)/p and |^ = ^ = O , or c 2 V4./l° and §£ ' ° ' 

The first case corresponds to longitudinal waves, in which the motion is 
along the direction of propagation and in the second case the motion is 
transverse and parallel to the wave front. 

The theory cf transverse elastic body waves was first investigated 
by Navier (1827) and, a little later, more rigorously by Poisson (1829). 
These treatments appeared about the same time as Fresnel's theory of the 
transverse nature of light vibrations. Since prior to this the concept of 
transverse vibrations propagated through a medium had not been con- 


sidered, subsequent developments in the theory of elastic waves tended 
to become associated with discussions on the propagation of light. 

(See for example Stokes [1848] and Kelvin [1904] .) 


The velocity of the shear waves depends only on the density and 
the shear modulus of the medium and it might appear intuitively that the 
velocity of the longitudinal waves should depend only on the density and 
the bulk modulus, k, but k = A^+ from (A. 8) and the shear modulus 

is also involved . The physical reason for this is that in the propagation 
of longitudinal waves the medium is not subjected to a simple compression 
but to a combination of compression and shear. To see this consider a 
small cube of material in the path of such a plane wave traveling in the | 
x-direction; its cross-sectional area normal to the x-axis will not alter 
during the passage of the wave but its x-dimension will be altered. Thus 
there is a change in shape of the element as well as a change in volume 
and the resistance of the medium to shear comes into play. 

Equations (A. 24) , (A. 28), and (A. 29) are all of the form: 


a 2 * 

at* 


c 7 v * 


(A. 35) 


and when the deformation is a function of only one coordinate , for example 
x, the equation becomes: 


lk$. 


d cC 2 

W = C 

The general solution for this equation is: 

qC = -fu-ct ) + P(X+ct) 


(A. 30) 


(A. 37) 


where f ana f are arbitrary functions depending on the initial conditions . 
F corresponds to a plane wave traveling along the negative x-axis and 
f to one along the positive x-axis. For each wave it may be seen that, 
if at any time t^ , c C is a given function of x at a later time , t 2 • it will 
have the same shape displaced along the x-axis a distance c(t 2 - t^) 

If the disturbance is spreading from a point, the deformation 

will depend only on the value of the radius vector from the point , ~t . 

2 2 2 2 
Since r = x + y + z we have: 


— 2 if . J_ 


bx‘ r ‘ dr 1 + 
With similar equations for 


( r z - x* \ ax 

l T z J dr. 


(A. 38) 


a«c 


and g so (A. 36) becomes: 


l 

d 


t* c l5r* ^ar/ 


a* 

a7* a*' 

c?(r<S-) _ _ 2 3 2 (r*ct) 


or 


at* 


- c 




(A. 39) 


This is of the same form as (A. 36) and its solution is: 
rcc «f(r-ct)+F(r*ct) 


(A. 40) 


S 


Now f and F represent spherical waves , f is diverging from the origin and 
F is converging . The amplitude is inversely proportional to the distance 
r in both cases. 




i 



S lk 9* 

PRECEDING PAGE^BLANK MOT FILMED. 

THE REFLECTION OF SOUND FROM A PLANE INTERFACE 


B-l. Reflection and transmission at a fluid-fluid interface 


Although no fluid-fluid interface exists in the system to be con- 
sidered in this report, this is the case which has both a, reflected wave 
and a single transmitted wave generated at the interface . By not having 
to deal with the generation of shear waves the problem is simplified con- 
siderably for a better intuitive understanding. The case of fluid-solid 
will then be covered to complete the system . The method used for both 
of these cases follows basically that of Brekhovskikh (I960) with some of 
the notations changed to be more consi stant with the rest of this report. 

The sound field will be characterized by the acoustic potential <f> . 
The particle velocity and the acoustic pressure in a harmonic wave (the 
time dependence is e + lu>t ) will be expressed through <f) by the equations 

v = - grad <f> , p = + icjp (fr (B.l) 


Consider the problem of the reflection of a plane sound wave at a 
plane boundary separating two media. The density of the medium in which 
the wave is incident , the upper medium , is denoted by p , and the 
acoustic velocity by c. The corresponding quantities in the lower medium I 
are denoted by p ^ and Cj. The angle of incidence is denoted by B and 
the angle of refraction by 6 The normal to the wave front lies in the 
plane of the diagram (Fig. B-l). 

4- i UJ t 

Using these notations and suppressing the time factor, e , 

v * 

the incident wave can be written as 


d 

/ inc 


= A exp [ik(-x sin Q + z cos 0 )] 


(B.2) 


where A is the amplitude of the wave . The reflected wave can be 
written as 


(i f = RAexp[ik(^x sin0-z cos© )] 


(B.3) 
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j s '--. v 
• **■ 

where R is the reflection coefficient. Assuming a linear system the -total 
potential in the upper medium will be 

+ ^ref s Aexp[ik{-xsin @ +zcos£)3 + RAexp[ik(-xsin Q +zcos©)] 

r: Aexp (-ikxs in Q) [exp (-f-ikzcos Q ) + Rexp(-ik 2 cos © 2 (B . 4) 
The refracted wave can be written in the form 

0 [ - TAexp[ikj(-x sin^ + z cos©j)3 (B.5) 

where T is the transmission coefficient and - ou /c^ is the wave 
number in the lower medium . 
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The boundary conditions for this system are that the acoustic 
pressure and the normal component of the particle velocity must be con- 
tinuous across the interface at z =0. Using (B.l) these conditions can 


be written as 




(B.6) 


a* 




(B.7) 


from which the two unknown coefficients R and T, as well as the angle of 
refraction ^ , can be found. r 

Substituting (B.4) and (B.5) into (B.6) gives 


p (l+R) = f, T«#p[-i(/(, *] 


(B.8) 


Since the left hand side is independent of x, the right hand side must also 
be independent of x for the equality to hold for all x. This yields the 
well known refraction law since the only way to satisfy the condition is 
for 


k sin e K Sin 9 , (B. 9 ) 

This relation expresses the equality of the phase velocities of waves in 
both media propagating along the interface. It can also be written as 

si* & a- u t c * 

5 j#0, ~ ^ where yi = (B.10) 


Applying this to (B.8) yields 

x . A 

"~P = Wl I +■ R ) 

Moreover, substituting (B.4) and (B.5) into (B.7) gives 
CCS© ('-«)*-■ " Tc “ 01 



(B . 12) 
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From (B.ll) and (B.12) comes 
m cos B - n cos B, 

I' ~ yn cos B + n cos Q t # 

Which can also be written as 

m cos fP “ T/fo*- 5/'n*0 


(B. 13) 


R • 


YY\ 


cos & + i/n*-s 7 h*© 


(B. 14) 


Thus the unknown quantities in (B.3) and (B.5) can be determined 
from the boundary conditions and the problem can be considered solved. 
Before progressing to another interface however, we shall consider some 
of the implications of these results . 


At normal incidence (0-0, - 0° ), (B. 3) gives 

p, c . - 0c 

p _ m - n _ r l 1 

m + n f>i c i + pc , 


i.15) 


The quantity z - pc is called the characteristic impedance of the medium. 
Using impedance the reflection coefficient can also be written as 

z i ~ z ? 


R = 


1 


+ z 


(B. 16) 


If the normal impedances and 2 are defined as 


P< C 


t i 


7 = 

■ -l £0 SQ, f 


z * -fis. 

C.S © 


the reflection coefficient for any incident angle becomes 


R = 


z i - z 

*1 ^ Z • 


(B. 17) 


As follows from (B.14) the reflection coefficient will become zero 
at an angle 0 satisfying 


wt cos 


Q = -yyi*- sin Q - 0 # 


(B. 18) 


In this case there is no reflected wave and the boundary is completely 
transparent. Solving (B.18) for this angle, Q gives 


s,r| e b = 


! p^ r 


(B. 19) 
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This angle is analogous to Brewster’s angle for electromagnetic waves . 
Note that 0^ need not be a real angle and complete transmission will not 
necessarily be observed in any particular case. Rather it is seen that 
the condition 


O* 




must be satisfied. Hence, when m >1, n must be such that 1 < n <m, and 
m <. 1 required 1 > n > m . 

When n <1 and the angle of incidence satisfied the condition 
sin 0 > n, total reflection will occur. In this case (B.14) becomes 

R , ™ IS E (B.20) 

yy\ CoS @ + i V 3 i TA & ~ VI * 


B-2. Reflection and Transmission at a Fluid-Solid Interface 

Now we shall generalize the problem considered above to include 
the case where one of the bounding media is a solid . 

The particle velocity at any point of a solid medium can be 
expressed through a scalar and a vector potential, using the equation 
(See Appendix A) . 

Tjj - - grdt) (j ■* curl v* , ,. f (B.21) 


In the special case of a plane problem oriented such that all quantities 
depend only on the coordinates x and z, and that the particle trajectories 


also lie in the xz-plane, the potential can be chosen such that only its 
y-component , which will be denoted by ^ , differs from zero . Then , 


according to (B . 21) , V will be a vector with components 

</” . 2d 7) 

” + i7 


a r* - 

K 


(B .22) 



15 ^. 


(j> and ^ are the potentials of the longitudinal and shear waves , 
respectively. In Appendix A these potentials were shown to satisfy the 
wave equations 

£ 

dt 1 (B.23) 

where 



c, * i0>» * 2 j», Vp 

(B. 24a) 


(B.24b) 


are the velocities of propagation of the longitudinal and shear waves as 
before . g 

The normal components of the stress and displacement must be 
conti nu ous across the boundary between the solid and the liquid . The 
tangential component of stress must also be continuous, but since the 
tangential stress cannot exist in the fluid this condition reduces to the 
requirement that the tangential stress at the boundary of the solid be 
ze**o. 

In the case of the plane problem the following stress components 
are of interest (A. 7) 

_ x ($j± + + p . ^ r\ 

^ 2 ." 3* 3* / ^ ^2 (B .25a) 

fh 

it 




where u and w are displacements along the x- and z-axes, respectively, 
as before. The z^axis is taken as normal to the interface. It is useful 
to express the displacement and the stress in terms of the potentials ^ 
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and . To do this use (B.22) and take into account that the displacement 
components u and w are obtained from the velocity components v and v 

X z 

by dividing by + i uj . Quantities referring to the solid will be denoted 
by the subscript 1, and the quantities referring to the liquid will have no 
subscripts. The elasticity of the liquid will be characterized by 
the acoustic velocity c, and the density f using (B.24a), with ^ set 
equal to zero . 

The sound field in the liquid is characterized by the potential <j> . 
Obviously, all the relations obtained for the solid medium can be extended 
to the liquid by setting - 0 and /-t = 0 . In particular, from {B.21) 
the connection between v and (f> will be of the form 

v = -grad (fa (B. 26) 

Then from (B.l) the acoustic pressure will be 

(B.27) 
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where , again , A is the wave amplitude and O is the incident angle . 

The reflected wave may be written in the form 

(B>32) 

Thus , the total sound field in the liquid will be 

j^Aje* expC-ik* s '" 6) (b 33) 


A longitudinal and a shear wave will be present in the solid. 


These waves can be written as 
rh - A Tetp\-iK, (t 5 

n L (B.34) 


if, = A 


(B.35) 


where k, k^ # andJ^j are wave numbers 


b r l< * J£- *.»< i 

'< C 9 I C, 1 


a? 


c. 


(B.36) 


and and Y ^ are the angles between the z-axis and the normals to 
the wave fronts of the longitudinal and shear waves in the solid. These 
angles are defined in Figure B-2 . 
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Figure B-2 . Geometry of the Fluid-Solid Interface 

The problem can be considered solved when R, T, P, 8 ^ and ^ j 
are all known. Substituting (B.33), (B.34) and (B . 35) into (B.28), (B.29) 
and (B.30) with z s 0 yields three equations from which these quantities 
may be found|. For example (B.30) gives 

k coS 8(3-0= “K 4 Co5 S/w kSi»&)*\ 


(3.37) 
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Since the left hand side of this equation is independent of x, the 
right hand side must also be independent of x. This can be true only if 
the equation 


k sin & * K Si* § Sin 


is satisfied, whence the directions of the waves in the solid are 
determined. 

Now (B.37) can be written 


K oos 


& (R-0 r “ k, <-oS &i~7~+ S i » R 


Similarly 


from (B.29). Adding and subtracting <1 aJ __ — on the right 

r ** ' n . 2. 


side of (B-28) and using 

ZL - tA 

™ e> X - 5 * 

this equation can be rewritten as 




Remembering that 2 

^ - pc 2 = f 9 ^ K 2 

and the wave equation 
(B.41) can be simplified to 

££) .t.-o 

m Y ™ fc. / 


2 

-7* 




(B.42) 
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where 


W " ~p (B.43) 

Substituting the values of $ , $ j, and into (B.42) yields a 

third equation for the determination of the coefficients R/ T and P , 


m ( , + ^) ~ (l-2 Tz siy' 2 &)~T ~ SivZf'P (B.44) 

A, 

Solving the system of equations (B.39), (B.40) and (B.44), and 
making some transformations using (B. 38) gives 

r _ z > t ZsIzlIIl: ~L 

t.s 2 Y, si v72 + ZL 

~ z,- si, J e[4z 5 - ezj - ± f +4Z j + z <B - 45) 

where 2, Zj,and Z g denote, respectively, the impedances for sound waves 
in the liquid and longitudinal and shear waves in the solid 


Z= 


_f e _ 

ccS& (7\ 


~7 = ft C l 
1> . cos (9 a 


Mof 7 

4- o Cd 9 i 


(B.46) 


_ I 2Z, c.i 2Y, 

W Z'Ca^ZY' *-Z s S)v72Y,*Z 

-± gz, 67-2^ s;»*e) 

_m Z,-f|s,Ve[4z r 5z,]- 

_._i 2Z S 5/w 

m z CO*? 2. y, Zj S » 11* y, 4-Z 

4Z.lr sinQ7fc*-c|v'"*e' 

2,-J* si/ef4z s - «z,]- ^ yi/sfz *42,7+* 


(B.47) 


(B.48) 
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Now, let us consider some implications of these results. At 
normal incidence (© » ^ ^ 1 ~ have ' 

Z, -Z — I. 2Z, j p = 0 



(B.49) 


which corresponds to the fluid-fluid case since no shear wave is generated. 
On the contrary setting T * 0 means 

- 2 “f 51 y \ 2 & “ 0 


or 




which gives 1 - 45° since 



from (B.28) we have 



hzl 

z s + z 


/ T - 0/ J P 


J- ^ 


(B.50) 


implying that shear waves but no longitudinal waves will be generated in 

... 'U \ 

the solid. 

Now consider the case in which a longitudinal wave of amplitude 
A is incident from the solid onto the solid-fluid interface. This wave will 
excite the following system of three waves at the boundary: 

(1) a reflected longitudinal wave (amplitude B^) 

(2) a reflected shear wave (amplitude B ) 

(3) a sound wave in the liquid (amplitude D) . 


l6l 


The entire system of waves may be written in the form 


(f) -O + 


for the sound wave in the liquid 


(() - A C X p [?^ < X sin &, + £ c o S 0 , jj 

•f B €Xp£i-k(x S/*i O.-ZCoS & f jJ 


for the incident and reflected longitudinal wc ves and 

% * & t exp[-i4 t (x si*r,-*eosY,jJ 

'•/’ 

for the reflected shear wave . 

Using the boundary conditions (B.28), (B.29) and (B.30) as before 


gives 


d'= Z*Z*lil 2 A ~ 

^ - Z * Z siSZl + Z, 




c* . 4, 


_ Z* ^5-" g.(<Z > ->-5-Z,)-4|i 5lVO. (Z, ■‘• Z, ) - Z, 
” z- ^S.- a 6>(4Z J -?Z,)-4^ s,^(2 s -Z,) + Z, 


-y'=(i’R') — . - c ^ 

c , cos 9 cos 2 8.X 


= 0-R') 


c cos 


c*s/»*o;jj 7 5S 5 /Ve, + s^7 


(B. 55) 
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5 I v» ZQ, 

cos 2t { 



c| 5 in 2 O, 
c*- 2 c* s ,V<S> . 


(B. 56) 


Similarly when a shear wave (of amplitude A .) with particle motion 

s 

in the xz-plane is incident on the boundary from the solid the system of 
waves can be written in the form ’ 


<p=D e*p[- i/<(xs/'»£>*a«ose)] 


(B.57) 


for the sound wave in the liquid ,, 

(y Q t -Z cos ^ 

for the reflected longitudinal wave in the solid; and 

0 : i 

- ^ = A $ e*p{~L-Ji t (y s / m ^ * 2 cos % j] 

vi + B s e y p iM t (x si * ^ - 2 c a s )J 

for the incident and reflected shear waves . 

Again using the boundary conditions yields 

= n u - z *^' *••*&}* z * *'•** 2 K 

■ ■ " z t B *Zi,+Z % Si» !! 2t l 


(B .58) 


(B. 59) 


‘v 


(B. 60) 


w 
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rV’N t v>© 

T " 0 + R > TT^X~ 

— (1 B ) 2 aiw % jfcf-c'sisi'jr' 


P - -£)V"") 


cos P ^ 

S/n 2e, 


- ~(i*R“) irH'ZA 

V 2 S/71 ^ 7c| - c* 


57^r 


(B .61) 


(B. 62) 


The only case re m aining to be investigated is a shear wave with 
particle motion parallel to the interface. A wave of this type creates 
neither normal displacements nor normal stresses at the interface. This 
means it will not excite a sound wave in the liquid or a longitudinal wave 
in the solid . The one boundary condition left is satisfied by a reflected 
shear wave of the same amplitude . 

■q Thus all waves incident on a fluid-solid interface have been 4 
considered since, in an isotropic, linear solid, a plane wave of arbitrary 
polarization can be considered as a superposition of waves of the polari- 
zations investigated. 
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EXPERIMENTAL DETERMINATION OF THE MATERIAL CONSTANTS 


C-l . Wax Floating on the Water 

The piece of wax used was cast with both sides flat and parallel. 
The thickness was then measured and found to be 37 mm. The transducers 
were aimed at the surface of the water from below and carefully adjusted 

jtr 

for maximum signal reflected from the smooth water to determine normal 
incidence . This signal was taken as coming from a perfect reflector and 
recorded for comparison with that reflected from the wax. The wax was 
then floated across the illuminated area and the return observed. Two 
distinct pulses were found and the amplitudes and relative timing were 
recorded. With this information it Should be possible to obtain an estimate 
of the reflection coefficient/ the attenuation constant, and the velocity 
of propagation for the wax . ... { 

The measured values were: 


Return 

Calibration 
First pulse 
Second pulse 


Attenuation 

70 dB 
46 dB v 
46 dB 


Voltage Time Delay 

3 vpp 203 l^s 

2.4 vpp 1994 

2.2 vpp 2032^8 


''Attenuation" refers to the value of attenuation inserted ahead of the 
receiver to avoid saturation and the difference corresponds to the change 
of system gain. Using these values, (B.49) and (B.55) the desired quanta 
ties were calculated as follows: 


v = 


-1 

At 


2( 3,7 x f O' 2 ) 
38 X I O’* 


190 0 r»/$ec 







(C.l) 


(C .2) 
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where V caJ is the signal received from a perfect reflector and V^ gt is the 
amplitude of the first pulse from the wax. Since the system gain was 
different for the two voltages as measured this must be taken into account. 
This gain change was 24 dB which corresponds to a voltage ratio of 1/16, 
yielding: - 

R = 3Ti t- 05 . (Ci3) 

Z -2 

From (B.49) R - z'+z. whicu for Zj> Z gives a positive reflection 
coefficient. Since Z » pc and Zj P\ c i this value can be calculated as 
a check by 

Z * 1.5 x 10 6 Zj =• (.9) (1900) x 10 3 = 1.7 xl0 6 

R = ~flT ~ 0,062 (c.4) 

which is in good agreement with the measured value . 

The attenuation constant, oC^ , was determined by using the 

relation 

V ! V,TTV ‘' 1 (c.s) 

<2^ caf. c 

where T = *f; (1-R) is the transmission coefficient from the water 
to the wax, 

T- — (1-R') is the transmission coefficient from the wax to 
C 1 the water, 

x is the total distance traveled through the wax, and 
V 2nd is am Plltude °* the second pulse . 

This is the equation for transmission through a section of material but 
since the wax/air interface has a reflection coefficient of unity it ?/ V 
applies to this case. Solving for and substituting the experimet^illy 
determined values yields 
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TTV 




vi 


cal. 


v,._< 


< / ( /.oo *U0 3)(3X'6) 

“ 7.4- /vt 2.2 


“ 7A ^ 1 ^ “ #33*<f/ t 


(C . 6) 




These results were such that it was thought the wax was worthy of further 
testing . 


C-2 . Longitudinal Waves Through a Section of Wax 

For a more accurate determination of and a rod of wax 
was cast with a diameter of 50 mm and a final length of 128 mm after the 
ends were cut square. This is thick enough that the velocity of propaga- 
tion is very close to that of an infinite piece . 

The transducers used were first placed in direct contact with each 
other, with a small amount of wax applied as a couplant to insure tight 
coupling, and the received voltage was recorded. They were then applied 
to the ends of the rod and the time delay and attenuation introduced by 
the rod were measured to give 

&t - 68^ s 
V o /V i = 1/312 

These data were used to calculate the following constants 


v A = 128/68 = 1870 m/sec 

(C. 7) 

In 312 - = .45 nep/cm 

(C.8) 


which agree very well with the preliminary results. 
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C-3 . Wax Prism 

To measure the shear wave velocity a prism was constructed of the 
wax. By illuminating one face of the prism at an angle other than normal 
incidence a shear wave as well as a longitudinal wave was generated in 
the wax. At the other side mode conversion would also occur, generating 
a longitudinal wave in the water at an angle determined by the refractive 
index between the velocity of propagation in the water and in the wax. 

The prism was cast with an included angle of 45°. The available 
molding facilities would only allow the prism to be 20 cm from apex to 
base, which restricted the angles that could be used without having a 
corner interfere. The faces of the finished prism were not as flat as 
desired due to the shrinking of the wax on cooling but the curvature was 
determined to be small enough that a fair estimate of the velocities 
could be made. .o 

This prism was submerged in the water with one transducer aimed 
at the center of the front face at 22 .5° from normal incidence (perpendicular 
to the center line of the prism) . The second transducer was then scanned 
in an arc behind the prism and the signal received was recorded as a 
function of angle. Two distinct peaks were found as expected. The longi- 
tudinal wave gave a large peak at 12° and the shear wave produced a 
much smaller peak at -38°. 

Referring to Figure C-l the equations for this prism are 


_ V 
n= 

(C.9a) 

sin ©j = n sin 

(C.9b) 

sin © 2 n sin 

(C,9c) 

it 

cn 

O 

(C.9d) 

e 2 = 22.S°-jZS m 

(C.9e) 


where n is the desired result and the angles 0^ and (j) are known. 


Figure C-l . Definition of Angles for Derivation 
of Prism Equations . 


These equations can be used to find n directly but the expression 
is very complicated so an intermediate step of solving for ^ was used. 
Start by substituting (C.9d) into (C.9c) which gives 

Sin Q t - n sin (45*- 

*n(sin45*c°*0 z - c 09 45* *i» ^ ) 
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r'i 


since sin 45°= cos 45° « (2) . Dividing both sides by sin yields 


S in 0, _ jn_ f cos 0?- sin <fi 2 
Sin Q z iZ \ sin © 2 


-) 


(C.ll) 


but using (C.9c) this can be written as 
s i n 0, l_/cos^ a 

S/v* 0 2 if? \ sin 0 t J 

= yf ( eot h ~ i) 


(C. 12) 


or 


r Z Sin 


(C.13) 


which can be used to find $2 s ^ce ^ j and ©2 were determined 
experimentally . Once i> 2 is known (C.9c) is used again tc find n by 


Y\ = 


si n &£ 


Si 


"A . 


(C . 14) 


Applying this to the experimental results of the longitudinal wave 


we get 


e 2 = 22. 5*-lZ'=l0.5‘ 

, r-r* sin 22. S' 

co t ^ - V? - w - Vo -^- *1 


3 & 3 

= 1.414— + | r 3.3Q 

.182 


(C. 15) 


Which means 


<j> • cot 3.98 ~ 14.1* 


(C* 16) 


giving 


. taz 

n = ^T^.745 
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c / So o , 

c ( = ? • 7 45 ~ - 2010 m/sec 


Similarly for the shear wave 


A * 2 2.5'-(-38') = GO. * 
__ s i v* 22 . , 

cot <t> - -JY 1 


= '-* ,+ .870 +l = ' 633 


^ = co 


i'7.G33 =3/.3' 


.07g _ 

.5-22 




V /. 66 


= 900 w/sec 


Thus an estimate of the velocities was made. The longitudinal , 
wave velocity found was higher than that previously determined but was 
of the same order of magnitude , with the error being attributed to the 
curvature of the prism sides. This error should also be present in the 
measurement of the shear velocity which would imply an actual velocity 
somewhat less than the 900 m/s determined here. 
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04 . Method of Hughes 

The method of Hughes (IP 49) consists of driving one end of a right 
circular cylinder with a pulse -modulated carrier of longitudinal waves 
and detecting longitudinal waves at the other end. The detected signal 
consists of a series of pulses with the delay times the important para- 
meter. The delay from the start of the transmitted burst to the leading 
edge of the pulse string , T , corresponds to the time for the longitudinal 
wave to traverse the material under test . The time between the first and 
second pulse. At, is the sum of infinitely many waves that have tra- 
versed a path such as ABCD shown in Figure C-2. These arise because 
the transmitted longitudinal wave is not strictly parallel to the surface 
and mode conversion takes place at the reflection. The path AB is 
traveled as a longitudinal wave, BC is traversed as a shear wave, and 
CD is again a longitudinal wave. All these waves suffer the same delay 

for crossing the rod one time in shear so the result is a second pulse at 

u 

' the detector end. Also present are following pulses generated by paths 
such as ABCEF. 


DRIVEN 

END 



Figure C-2. Two Possible Paths for Sound Waves 
in a Solid Cylinder 


As shown by Hughes the time between pulses is 

% 


At =Dp?-<? J. 


(C . 24) 


,// 


This equation can be solved for c 2 since the diameter of the rod, D, and 
o, are known quantities and A t is measured experimentally. This 
gives 




This method was implemented by casting a circular cylinder of 
the wax then mounting a transducer on each end of it . This first attempt 
did not produce the distinct pulses as expected but only an extended 
trailing edge of the main pulse . The cylinder was then cut open to 
determine the cause of the discrepancy. It was found to have large 
cracks at the center caused by the stresses produced by contraction of 
the wax as it cooled in the mold. 

To eliminate these cracks it was decided to extrude a cylinder for 
the test . A mold was prepared with one end Slightly smaller than the 
main section. This was filled with molten wax then slowly cooled while 
a reservoir of molten wax was maintained on top to reduce the formation 
of stress cracks . After being thoroughly cooled the cylinder was forced 
out the smaller end to form the final test piece. The ends of this were 
then sectioned and determined to be of a homogeneous nature . The trans- 
ducers were then mounted and one was driven with bursts of 350 KHz with 
pulse length of 60 ^us . The signal received at the other end was then 
observed and the pertinent times recorded . Photographs of the received 
signal are shown in Figure C-3 . //_ 


I'O : 

t:' !' 


, I! 


i-’ 


176 . 



(a) Start of Trace Corresponds to Start of Transmitted 
Pulse. Time Scale = 50//is/cm 




(b) Expanded Sweep Showing At More Clearly. 

Time scale » 1(3/4 s/cm. 

Figure C-3. Photographs of Signal Received at Detector 
End of Rod Number One. 
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A second rod with a different diameter and length was then con- 
structed by the same technique to serve as a check of the experimental 
method. The different dimensions would eliminate any possibility of a 
spurious spatial resonance effect which could have caused erroneous 
results . The dimensions and experimental results for both rods are 
tabulated below: 


Parameter 

Rod #1 

Rod #2 

Length 

22.5 cm 

13.0 cm 

Diameter 

3 .9 cm 

4 4 cm 

Frequency 

350 KHz 

350 KHz 

T 

125>tS 

73 /js 

At 

68 /as 

76 /as 


1800 m/sec 

1782 m/sec 

C 2 

550 m/sec 

552 m/sec 


The good agreement between the two independent experiments should not 
be considered as an indication of true accuracy but it does give more 
confidence in the results obtained , 


G-5 . Shear Attenuation 

To measure the shear attenuation it was necessary to excite a 
shear wave through a section of the material and determine the loss per 
unit distance . This would have been easy if shear mode transducers had 
been available. Sui£ft was not the case, however. Shear waves are also 
generated by mode conversion when an interface is illuminated at other 
than normal incidence, as shown in Appendix B. To eliminate errors due 
to inaccurate knowledge of the material parameters, it was decided to use 
a difference method by comparing the change of attenuation between two 
different thicknesses of the material. To reduce the effect of the longitu- 
dinal wave , the interface should be illuminated at an angle beyond the 
critical angle for these waves . The critical angle can be determined by 
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setting the angle of refraction equal to 90° and using Snell's law to solve 
for the incident angle. This gives 


€j . = sin (f; si nS, ) 

f.r 0 , = 9o* 

-i /r \ . -»/ 1 5*00 \ 

Or s '" (cj s s,n \T&o°) 


- 5 6.5*: 


(C.26) 


Using an incident angle of 66° should thus give no longitudinal waves in 
the material . 

Two slabs of wax were prepared with thicknesses of 1.4 cm and 
4.3 cm . These were mounted end-to-end such that the transducers 
could be moved from one to the other without changing anything except 
the thickness of wax. The difference between the two thicknesses was 
53 dB attenuation as measured experimentally. 

The path length of the shear wave through the wax can be deter- 
mined from the incident angle, 0 , Snell's law, and the thickness, T, by 

= S/» (^SiyiQ) 


(C . 2 7) 




T 


T 


**• ~ VT-fisiuaf . 947 , 


For the thinner piece 


L = 


- 1.4 


.947 


- 1 .48 cm 


and the thick section 


.W * 4 - 54 cm . 


Thus A L - 4.54 - 1. 48 =3.06 cm 
and the attenuation per centimeter is 


Otfii ' 


At t VI ^ - 

cm 


TC - TM ' 17 - 2S dB/om 




(C .28) 


(C.29) 


(C.30) 


(C .31) 
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for the shear wave in wax. This corresponds to 


oC - 2.75 nep/cm. 


C-6. Material Constants 

With the values of c, , C^and p determined by the preceding 
methods the results derived in Appendix A can be used to determine the 
mechanical constants for the material. Lame's constant was calcu- 
lated by using the relation 

which means 


^ , pci : (9.37)Cf.tro)*KI0 < ' 


- 2 .83 x 10 8 newton/m^ . 


Similarly A m can be found from^,-^ f p and g, by 


c , * 




Solving for A m gives 


=^9,37 K'.8) t -2(289)]l°* s 2.4-7 X I0 S ( C< 


Once Lame's constants have been determined any of the other 
commonly used elastic constants can be calculated from (A. 6) to (A. 8) . 
young's modulus, E is found by (A. 6) to be 

r s „ a o i j. / m 

C. — — * r»«wr*w/wiV. . 


i 
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Poisson's ratio is 
hy* 


-i = 




= . 449 


from (A. 7) . Finally from (A. 8) the bulk modulus k m is 
^m ” + 3~ ~ 2.65 * 1 0 9 new t oy s 

C-7 Wax- Water Interface 


m 


(C . 3 5) 


(C.36) 


The equations of Appendix B were programmed to allow the computa- 
tion of reflection and transmission coefficients for the three possible 
cases, i.e. longitudinal wave in water hitting the wax, longitudinal wave 
coming out of the wax, and a shear wave striking the water from the wax. 
Also computed were the magnitudes of the two-way transmission through 
the wax. This program computed these quantities for each case with the 
angle of incidence changed in increments of 5°. The only ones found to 
be significant were the longitudinal waves. This program is reproduced 
in Figure C-4. The coefficients of interest appear in Chapter III as 
Figures 3-11 — 3-16. 
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ARSIN(X)--ATAN(X)/SQRT(1 .o-x*x) 

C =1500 

Cl=1800 

C2-550 

RHO=1000 

RHOl- 937 

FM= RHOl/RHO 

RC-RHO*C 

R1C1=RH01*C1 

R1C2 = r RHOl*C2 

ClOCl/C 

C2C=C2/C 

C1C2-C1/C2 

C1C22=C1C2*C1C2 

WRITE (6, 5) 

5 FORMAT (3X , 5HTHETA , 8X , 1HR , 1 OX , 1HT , 1 OX , 1 HP , 9X , 2HR1 , 9X , 2H 

1TI ,9X, 2HP1 , 9X,2HR2 , 9X, 2HT2 , 9X, 2HP2 , 8X, 4HT*T1 , 7X, 4HP*T2) 
DO 20 I 2,11 
THETA-(I-1)*5 
TH- THETA*Q .01745 
STH=SIN(TH) 

CTH-COS(TH) 

TTH--STH/CTH 
C1CSTH"C1C*STH 
C2CSTH «C2C*STH 
TH1» ARSIN(CICSTH) 

CTHl=COS(THl) 

TTH1=2.*TH1 
STTH 1“ SIN (TTH 1 ) 

GM1'ARSIN(C2CSTH) 

CGM1*C0S(GM1) 

TGM1 = 2.*GM1 
STGM1= SIN(TGMl) 

CTGM1»C0S(TGM1) I 

STGM12=STGM1*STGM1 ' 

CTG M 1 2 =CTG M 1 *CTG M 1 
Z*RC/CTH 
ZI = R1C1/CTH1 
ZS=R1C2/CGM1 
ZlCS t Zl*CTGMl 
Z1CS2»Z1CS*CTGM1 
ZSSN- ZS*STGM1 
ZSSN2^ZSSN*STGM1 
D : Z+Z1CS2+ZSSN2 
R t (~Z+Z 1CS2+ZSSN2)/D 
FMD-FM*D 
T-2 . *Z1CS/FMD 
P--2 . *ZSSN/FMD 
Rl = (Z-Z1CS2+ZSSN2)/D 
T1 = CTH1*(1 .-RI)/(C1C*CTH*CTGM12) 

P1 = STTH1*(1 .-R1)/(C1C22*CTGM1) 

R2 s (-Z-Z1CS2+ZSSN2)/D 
T2 = TTH* (1 . +R2)/(2 , *STGM 1 2) 


% 

w 
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P2- - (CXC22*CTGM 1 * (1 . +R2))/STTH1 
TT1=T*T1 

PT2?P*T2 % 

WRITE (6, 10) THETA, R,T,P,R1 ,T1 /PI ,R2 ,T2 ,P2 ,TT1 ,PT2 
10 FORMAT<12(1X,1PE10.2)) 

20 CONTINUE Y\ 

STOP \\ 

END V\ 


Figure C-4 . Program to Determine Reflection and Transmission 
Coefficients for Wax-Water Interface, 





NASA — MIC 


